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SECTION  I 
INTRODUCTION 


The  response  of  materials  to  impulsive  loading  is  a  criti¬ 
cal  area  of  research  in  this  country.  A  recent  overview  of  this 
subject  prepared  by  the  National  Materials  Advisory  Board  (NMAB) 
has  identified  the  limitations  of  materials  models  as  the  "prin¬ 
cipal  limitation  in  the  use  of  numerical  methods  that  are  essen¬ 
tial  in  the  design  of  many  ballistic  systems".  2UKAS  states  that 
"the  description  of  dynamic  material  behavior,  especially  mate¬ 
rial  failure  under  high  loading  rate,  remains  the  greatest  single 
obstacle  to  the  accuracy  and  utility  of  computer  codes  for  solid- 
solid  impacts."  The  NMAB  study  also  recognized  that  the  charac¬ 
terization  of  fracture  and  dynamic  plasticity  are  the  areas  most 
in  need  of  investigation. 

The  work  reported  here  concerns  dynamic  fracture  in  metals. 
The  fractures  occur  under  impact  and  are  of  a  type  commonly 
referred  to  as  spall.  In  this  introduction,  we  present  the  basic 
phenomenology  and  provide  a  theoretical  framework  for  analysis. 

The  experiment  design  is  discussed  in  Section  II.  The  subsequent 
sections  present  results  for  the  individual  metals  investigated  — 
copper,  mild  steel,  4340  steel,  armor  steel,  nickel,  and  aluminum. 

1.  IMPACT-INDUCED  TENSILE  FRACTURE 

High-velocity  impact  events  generate  intense  compressive 
shock  waves  that  propagate  away  from  the  impact  site.  Wherever 
the  shock  waves  encounter  free  surfaces,  release  waves  are  gen¬ 
erated.  Release  waves  communicate  the  f ree-surf ace  zero-stress 
boundary  condition  to  the  shock-compressed  material.  Interaction 
of  release  waves  results  in  tensile  stresses  that  can  cause 
fracture . 

The  nature  of  shock-induced  fractures  is  often  surprising 
and  counter-intuitive.  Some  examples  are  shown  in  Figures  1 
and  2  . 
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Figure  1. 


Examples  of  Fracture  in  Steel  Targets  Struck  by  11  mm 
Diameter  Balls  at  3  km/s  (from  BLESS  and  BARBER,  1979) 


Examples:  of  Spall  from  Present  Program.  Steel  Targets 
Struck  by  2  mm  Flyer  Plates  at  300  to  400  m/s.  (Spall  f 
a  flat  plate  (a);  and  spall  in  a  conn  struck  on  its  r  a 
(b) ,  from  shot  276.) 


Figure  1,  taken  from  BLESS  and  BARBER,  1979,  shows  cross 
sections  through  mild  steel  targets  struck  by  ball  bearings  at 
3  km/s .  It  can  be  seen  that  the  shape  of  the  target  free  sur¬ 
face  apparently  had  a  very  dramatic  effect  on  the  geometry  of 
the  fracture  pattern.  Radial  cracks  occurred  in  the  hemisphere, 
but  a  crack  parallel  to  the  impact  surface  occurred  in  the 
cylinder. 

Figure  2  shows  examples  taken  from  the  present  work.  Part 
(a)  shows  fracture  patterns  in  a  cone  struck  by  a  flat  flyer 
plate.  At  slightly  higher  velocities,  the  fractures  coalesced, 
so  that  material  separation  occurred  on  four  separate  fracture 
systems  --  two  axially  symmetric  and  two  radial.  Part  (b)  shows 
1-D  spall  in  steel;  a  flat  flyer  was  launched  against  a  flat 
target,  resulting  in  a  plane  fracture  surface. 

Spall  plays  an  important  role  in  many  processes  of  practical 
importance.  Two  examples  from  different  fields  are  shown  in 
Figure  3.  Part  (a)  shows  a  cross  section  through  a  piece  of 
titanium  struck  by  a  glass  bead  at  7  km/s.  If  not  for  the  spall, 
penetration  would  have  been  incomplete.  This  particular  specimen 
and  projectile  represented  a  model  of  a  meteoroid  impact  on  a 
spacecraft  (from  BLESS  and  GREEN,  1980).  Part  (b)  shows  a  small 
scale  model  of  naval  ship  armor.  A  projectile  has  almost  pene¬ 
trated  the  armor,  and  can  be  seen  protruding  through  the  rear 
surface.  A  spall  plate  has  been  detached  that,  in  a  full  scale 
impact,  would  measure  about  10  cm  across.  Such  a  plate  would 
be  quite  a  hazard  to  interior  components  of  a  ship. 

2 .  CONCEPTUAL  FRAMEWORK 

Ultimate  fracture  in  metals  has  been  a  distinct  field  of 
specialization  for  many  years,  and  a  great  many  results  have 
accumulated.  However,  very  little  of  the  literature  on  this 
subject  is  applicable  to  the  dynamic  bulk  tensile  fractures 
associated  with  spall  formation. 
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Figure  3.  Examples  of  Spall  Associated  with  Penetration. 

(a)  Cross  Section  Through  a  Titanium  Plate  Struck 
at  45°  by  a  1.4  mm  Glass  Bead  at  6.8  km/s  (from 
BLESS  and  GREEN,  1980) .  (b)  Rear  Surface  of  25  mm 

Aluminum  Plate  Struck  by  a  4  mm  Steel  Rod. 
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It  has  been  established  that  spall  fracture  is  influenced 
by  (1)  maximum  tensile  stress,  (2)  deviatoric  stress,  and  (3) 
stress  history.  The  only  current  model  that  takes  all  of  those 
factors  into  account  is  the  NAG  theory.  This  model  has  been 
described  in  a  series  of  reports  from  SRI  International  (see, 
for  example,  SEAMAN  et  al.).  There  are  several  disadvantages  of 
the  NAG  theory  that  have  discouraged  its  use  in  the  present  work. 
First,  being  micromechanical,  it  contains  a  very  large  number  of 
parameters  which  are  evaluated  from  relatively  tedious  experiments 
near  the  spall  threshold.  Second,  since  the  model  is  "calibrated" 
near  the  spall  threshold,  it  may  be  less  accurate  where  threshold 
criteria  are  greatly  exceeded  (the  usual  case  in  ballistic  impact.) 
Third,  the  NAG  model  is  only  incorporated  into  a  few  finite  dif¬ 
ference  codes.  In  particular,  at  this  time,  it  is  not  an  option 
to  the  sponsors  of  the  present  program. 


The  approach  here  is  based  on  separate  time-dependent  and 
stress-dependent  criteria.  Time  dependence  is  assumed  to  be 
described  by  the  TULER-BUTCHER  criterion. 


Kx  =  f{  |o-cry  j)  Adt 
K ^  >  K ^  for  spall 


(1) 


According  to  this  model,  there  is  a  minimum  normal  stress,  o^, 

at  which  spall  can  occur*.  At  o  . _  =  an ,  the  tensile  stress 
r  mi  n  u 

must  be  applied  for  an  infinitely  long  time  for  a  spall  plane 

to  form.  When  n  ,  <  a.,  spall  takes  place  at  shorter  times. 

min  u 

For  very  negative  values  of  o  ^  ,  spall  takes  place  very  rapidly, 
and  the  time  dependence  is  less  manifest. 

BLESS  and  BARBER,  1979,  discuss  the  fact  that  \  =  2  appears 
to  be  appropriate  for  steel.  Equation  (1)  also  has  an  attractive 
physical  interpretation  for  \  =  2.  Therefore,  the  value  of  \  =  2 
has  been  favored  for  analysis  of  data  in  the  work  reported  here. 

*The  nomenclature  used  here  is  that  stress  is  taken  as 
positive  in  compression.  However,  comparative  terms  describing 
tensile  stresses  always  refer  to  absolute  values. 
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In  practice,  it  often  turns  out  that  spall  can  be  described 

by  a  single  spall  stress  value,  a  .  Spall  occurs  when  o  is 

s  s 

exceeded.  J5UKAF  (1979)  and  BERTOLF,  et  al.,  discuss  the  utility 

of  this  approximation.  The  success  of  this  model  may  result 

from  the  fact  that,  in  many  impact  situations,  9o/9t  is  very  hiqh 

on  the  spall  plane  so  target  behavior  is  not  very  sensitive  to 

the  exact  value  of  a  . 

s 

Most  research  on  fracture  initiation  has  been  motivated  by 

structural  applications.  Rate  effects  are  seldom  considered. 

The  concepts  that  have  evolved,  however,  can  be  very  helpful  in 

characterizing  dynamic  fracture.  It  is  widely  accepted  that 

fracture  criteria  can  be  expressed  in  terms  of  the  plastic 

strain,  i:  ,  and  the  ratio  of  mean  stress  to  llow  stress,  o/Y, 

P 

See,  for  example,  HANCOCK,  et  al.  Figure  4a  shows  a  graph  of 
failure  data  for  steel  plotted  with  these  parameters. 

The  stress  trajectory  in  o  and  coordinates  is  a  function 
of  test  geometry.  In  a  conventional  uniaxial  stress  test, 


o  =  o/3 
x 

Y  =  0 

x 

o/Y  =1/3 


(2) 

(3) 

(4) 


For  a  hydrostatic  stress  state,  r  =0.  On  a  uniaxial  strain 
path,  using  the  approximation K  ,  =  constant,  u  =  constant,  and 

b 

Y  =  constant, 

c  ,  =  o/M  (5) 

X 

c  -  2r  / 3  (6) 

p  x 

o/Y  =  3Kf.  /2 Y  +  K/2m  (7) 

P 

Below  the  elastic  limit,  a  =  o,,-.  »  and 

max  H 

o/Y  =  K/2y  (8) 


Figure  4b  shows  graphs  of  o/Y  vs 
representative  of  steel.  OERA  is  the 


Cp  trajectories  for  values 
case  of  uniaxial  strain: 
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Compressor) 


(b) 


Figure  4 


Tension 


^r3cture  of  Steel  in  Terms  of  a/y  vs  e 
(a)  Fracture  Criteria.  Curve  is  Afte?' 
OSBORN,  1978.  Other  date  are  from  BLESS, 
et  al,  1979.  (b)  Stress  Paths  for  Steel 


compression  followed  by  release.  RB  shows  a  tensile  path  more 
hydrostatic  than  1-D  strain,  namely,  e2  =  ~  O.le.^.  RC  is  a 

path  more  compressive,  namely  e2  =  ~  SD  1S  the  path 

for  a  uniaxial  stress  test.  DCAB  is  the  hypothetical  fracture 
surface.  Uniaxial  strain  and  uniaxial  stress  tests  give  points 
A  and  D  on  the  fracture  surface,  which  are  rather  far  apart  (in 
reality,  much  further  than  shown  in  Figure  4b) .  The  stress 
states  most  representative  of  impact  events  have  not  been  well 
specified.  It  is  unlikely  that  they  will  be  straight  lines  on 
Figure  4.  Actual  impacts  usually  represent  diverging  stress 
waves.  Hence,  compression  is  less  hydrostatic  than  equation  (7) . 
However,  release  states  may  be  either  more  or  less  hydrostatic 
than  RA,  depending  on  the  curvatures  of  the  release  wavefronts. 

3 .  OBJECTIVES 

The  objectives  of  the  proposed  program  are  to  quantitatively 
characterize  spall  behavior  in  a  number  of  materials.  The  mate¬ 
rials  investigated  were 

• OFHC  copper 
•1020  steel 
•4340  steel 
•armor  steel  (RHA) 

•nickel  200 
•1100  aluminum 

4 .  APPROACH 

The  experiments  performed  in  this  program  were  mainly  1-D 
impacts.  Some  2-D  experiments  were  performed  on  steel  targets. 
The  principal  diagnostic  technique  was  a  velocity  interferometer 
P'TSAR) .  This  device  provides  a  record  of  target  free  surface 
velocity  history.  The  principles  of  VISAR  operation  are  de¬ 
scribed  by  BARKER  and  HOLLENBACH.  VISAR  analysis  techniques 
used  here  are  described  in  Appendix  C.  Impacts  were  also  ana¬ 
lyzed  by  the  SWAP  method  of  characteristics  computer  code.  This 
code  was  developed  at  Sandia  Corporation  and  is  described  by 
BARKER  and  YOUNG. 
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The  general  approach  employed  in  the  program  is  presented 

in  Figure  5.  The  steps  in  this  approach  are  explained  below. 

1.  Review  all  pertinent  results  available  in  the  literature. 

Data  should  include  spall,  uniaxial  tensile  strength,  strain 
rate  effects,  and  fracture  mechanics. 

2.  The  spall  threshold,  uc,  is  defined  as  the  minimum  im¬ 
pact  velocity  that  results  in  partial  separation  of  a  spall 
plate  in  a  1-D  test.  Spall  thresholds  should  be  determined 
for  at  least  two  geometries,  comprising  differences  in  flyer 
or  target  thicknesses.  These  data  will  permit  assessment  of 
stress  duration  effects. 

3.  Check  that  there  is  indeed  a  well-defined  spall  threshold. 

In  some  ductile  materials,  considerable  damage  may  accumulate 
without  separation  of  a  spall  pla’  . .  In  some  others,  onlv 
partial  spall  plate  separation  may  occur.  If  a  material 
is  mainly  considered  for  armor  applications,  probat ly  spall 
plate  detachment  is  the  most  relevant  criterion.  However, 
in  a  charge  liner,  extensive  void  formation  is  probably  a 
useful  criterion. 

4.  In  the  event  that  a  spall  criterion  is  not  self-evident  or 
if  complete  separation  is  deemed  to  be  an  inadequate  measure 
of  spall  damage,  an  alternate  approach  must  be  developed. 
Usually  the  most  useful  approach  is  to  estimate  the  void 
content  of  the  spall  plane.  Void  formation  and  growth  may 
turn  out  to  be  a  threshold  phenomenon  that  could  be  used. 

The  effect  of  voids  on  strength  and  elastic  properties  may 
also  suggest  a  critical  void  concentration  that  cannot  be 
tolerated. 

5.  The  one-dim  visional  SWAP  method  of  characteristics  code  must 
be  calibrated  to  reproduce  the  observed  VISAR  records.  The 
best  literature  estimates  for  constitutive  properties  should 
be  used.  Of  particular  importance  are  yield  and  release 
properties.  If  these  are  not  available,  are  not  adequate,  or 


Figure  5.  Flow  Chart  Showing  Approach  for  Spall  Characterization 
of  Materials. 
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turn  out  to  be  inaccurate,  VISAR  data  unaffected  by  spall 
can  be  used  to  develop  improved  constitutive  properties  for 
the  undamaged  material. 

6.  Spall  should  be  introduced  into  the  SWAP  code  in  order  to 
reproduce  the  observed  spall  plane  and  the  observed  time  of 
arrival  of  the  spall  signal.  The  results  for  amplitude  of 
the  spall  signal  must  be  checked.  If  agreement  is  poor, 
then  the  material  properties  must  be  affected  by  incipient 
damage  in  a  way  not  accounted  for  by  the  SWAP  constitutive 
relations.  This  information  must  be  incorporated  in  the 
fracture  model. 

7.  Estimate  of  the  spall  stress,  a  ,  and  the  residual  stress  on 

s 

the  spall  plane,  a^,  directly  from  the  VISAR  free  surface 
velocity  records. 

8.  Develop  a  model  describing  spall  threshold  condition  based 

on  threshold  data,  VISAR  analysis,  and  the  SWAP  code  analysis. 
SWAP  will  allow  accurate  calculation  of  the  K^c  parameter. 

If  equation  (1)  is  inadequate,  then  develop  an  improved 
version. 

9.  Carry  out  experiments  well  above  the  spall  threshold.  Obtain 
VISAR  data  for  these  impacts. 

10.  Carry  out  SWAP  calculations  and  iterate  to  obtain  agreement 
with  VISAR  records.  Evaluate  spall  plane  conditions  leading 
to  rupture. 

11.  Use  VISAR  records  to  directly  obtain  estimates  of  ag,  a 
and  the  time  required  for  the  spall  p]ane  to  coalesce. 

12.  Determine  whether  or  not  the  over-driven  spall  criterion  is 

same  as  the  threshold  condition.  If  not,  introduce  a 

dependence  of  K^c  on  stress  rate.  Express  results  in  terms 

of  e  and  o/Y  also. 

P 

For  the  case  of  1020  steel,  step  12  was  followed  by  an  additional 
process.  Two-dimensional  targets  were  impacted,  as  shown  in 
Figure  2a.  The  complex  fracture  surfaces  that  resulted  were  com¬ 
pared  with  finite  difference  calculations. 
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SECTION  II 

SPALL  IN  1-D  IMPACT  EXPERIMENTS 


The  simplest  geometry  for  study  of  spall  phenomena  is  that 
of  "one -dimensional  impact."  This  term  refers  to  strain  in  the 
material  under  investigation: 

=  AV/V  (9) 

e  =  e  =0  (10) 

y  z 

One-dimensional  strain  states  are  commonly  produced  by  planar 
impact.  Many  discussions  of  planar  impact  loading  are  avail¬ 
able  (for  example,  ASAY  and  LIPKIN,  or  McQUEEN  et  al.).  The 
techniques  used  in  the  present  program  are  mainly  described 
in  Appendix  A. 

1.  WAVE  INTERACTIONS  IN  TARGET 

The  experimental  configuration  needed  to  produce  the  one¬ 
dimensional  state  is  shown  in  Figure  6.  A  flat  flyer  plate 
impinges  on  a  flat  target  with  negligible  tilt.  Observations 
are  confined  to  the  region  of  time  and  space  unaffected  by 
release  waves  from  the  edges  of  the  flyer  plate. 

Figure  7  is  an  x-t  diagram  showing  a  somewhat  simplified 
version  of  the  principal  shocks  and  release  waves  produced  by 
an  impact.  A  flyer  of  thickness  d  and  velocity  uq  strikes  a 
target  of  thickness  T  o.  like  material.  Impact  occurs  at  the 
origj-  .  Elastic  and  plastic  shock  waves  are  generated  at  im¬ 
pact.  The  elastic  wave  travels  at  the  compressional  sound 
velocity,  C  .  The  plastic  shock  waves  propagates  at  the  shock 
velocity,  which  is  nearly  equal  to  the  bulk  sound  speed.  All 
waves  in  Figure  7  are  represented  as  characteristics,  across 
which  the  Hugoniot  equations  apply: 

f-7_(U-u2)  =  p,  (U-u^  (11) 

*  °1  ~  P1U  *U2  ”  d2) 
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for  a  wave  travelling  into  state  1,  behind  which  the  material 
is  in  state  2. 

When  a  wave  carrying  a  particle  velocity  jump  Au  impinges 
on  a  free  surface,  the  surface  is  launched  with  a  velocity  in¬ 
crement  given  by: 

Av  =  2Au  (13) 

Shock  waves  reflect  from  free  surfaces  as  tensile  waves. 
Such  tensile  waves  are  termed  release  or  rarefaction  waves. 

In  this  example,  the  release  waves  are  also  assumed  to  behave 
in  an  elastic-plastic  way,  so  the  first  release  waves  propagate 
at  an  elastic  wavespeed.  Actual  release  waves  are  continuous 
because  release  waves  are  dispersive.  In  Figure  7,  as  in  the 
SWAP  code,  the  continuous  release  fan  is  modelled  as  a  sequence 
of  discreet  waves.  The  release  fan  from  the  flyer  plate  tra¬ 
verses  the  target.  If  the  target  is  relatively  thin,  it  may 
even  overtake  the  plastic  shock  wave. 

The  train  of  right-travelling  waves  arrives  at  the  target 
free  surface.  E  represents  the  elastic  wave  arrival,  P  the 
plastic  wave  arrival,  R  the  first  release  arrival.  To  the 
extent  that  the  first  release  wave  is  elastic,  the  time  of  R 
can  be  used  to  calculate  the  longitudinal  elastic  modulus  in 
the  shock-compressed  material. 

As  the  right-travelling  and  left-travelling  release  waves 
cross  each  other,  a  region  of  the  target  experiences  tensile 
stress.  If  the  tensile  strength  of  the  material  is  exceeded, 
someplace  a  spall  plane  develops.  This  is  labelled  RP  in 
Figure  7. 

2.  SIMPLIFIED  INTERPRETATION  OF  i RFE  SURFACE  VELOCITIES 

When  the  spall  plane  separates,  the  stress  relaxes  to  zero. 
Hence,  a  shock  wave  is  generated  that  raises  the  stress  from  the 
spall  stress  to  7cro.  In  general,  at  the  instant  of  spall 
failure,  release  waves  are  propagating  away  from  SP  behind  which 
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=  a  .  The  spall  shock  will  usually  overtake  a  portion  of 
s 

these  release  waves.  The  arrival  at  S  communicates  to  the  free 

surface  that  the  stress  in  the  spall  region  was  suddenly  reduced 

to  zero.  However,  since  the  shock  may  have  overtaken  several 

release  waves  on  the  way  from  SP  to  S,  the  magnitude  of  the 

shock  arriving  at  S  is  somewhat  less  than  a  .  The  spall  plane 

s 

presents  a  free  surface  to  incoming  waves;  hence,  later  release 
waves  reflect  as  shock  waves,  and  incoming  shock  waves  reflect 
as  tensile  waves.  The  free  surface  velocity  history  correspon¬ 
ding  to  Figure  7  is  shown  in  Figure  8.  Several  important  obser¬ 
vations  follow  careful  consideration  of  Figure  7  and  other 
variants  on  the  same  level  of  complexity. 


Figure  8.  Free  Surface  Velocity  History  Corresponding  to 
Figure  7. 
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It  can  be  seen  by  inspection  from  Figure  7  that  spall  takes 

place  first  on  the  tail  characteristic  in  the  release  fan  from 

P.  This  is  a  consequence  of  the  fact  that  the  two  intersecting 

release  fans  are  similar;  if  the  release  fan  from  the  flyer 

surface  was  much  more  compact  than  that  from  the  target  free 

surface,  then  spall  might  take  place  within  the  release  fan 

from  P.  Let  J  be  the  Riemann  invariant  along  the  characteris- 
P 

tic  leading  from  P  to  SP,  /» SP 

J  =  u /  -ttt-  «=  J  (14) 


Since  at  P, 


u  =  v 


then 


Hence 


-  v 

P  P 


Where  u  is  the  value  of  u  just  prior  to  sep; ration,  and  j 
s  p  L  s 

is  the  value  of  stress  just  prior  to  separation.  We  make  the 
assumption  p  =  constant  =  p  and  U  =  Cn ,  hence : 

O  o 

(’_  =  P  Cn  (u  -  v  )  , ,  r , 

s  o  B  '  sp  p  (16) 

The  J+  characteristic  at  SP,  the  instant  before  f:V  ure, 


[ 

J  =  u  +  I 

+  J 


Its  value  may  be  found  from  the  value  of  the  free  surface  velo¬ 
city  just  before  the  arrival  of  the  spall  plane  shock,  v  .  . 


J  =  v  . 
s  mm 
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Using  the  same 

assumptions  as  before 

o  =  p  C„  (v  .  -  u  ) 

s  o  B  min  3p 

(19) 

u  =  1/2  (v  +  v  .  ) 

sp  p  min 

(20) 

Finally , 

o  = 
s 

l/2p  Cn  (v  -  v  .  )  =  l/2p  C-Av^ 

'  o  B  p  min  o  B  s 

(21) 

A  derivation  somewhat  similar  to  the  present  was  given  by  DAVID 
et  al.  Equation  (21)  was  used  by  CHRISTMAN  and  ISBELL  to  deduce 
spall  stress  values  from  their  data.  According  to  them,  it  was 
first  proposed  as  an  empirical  formula  by  TAYLOR.  In  many  im¬ 
pacts,  p  c  may  vary  by  as  much  as  10%  between  the  target  free 
surface  and  spall  plane.  Therefore,  at  worst,  equation  (21) 
should  bo  accurate  to  about  *  5%. 

There  are  at  least  three  second-order  effects  that  modify 
the  above  portrayal.  The  first  is  associated  with  reflections 
from  discontinuities.  Wherever  two  characteristics  cross  each 
other,  a  discontinuity  results  across  which  there  is  a  density 
contrast.  Later  characteristics  incident  on  the  discontinuity 
may  spli ‘  or  reflect.  For  example,  when  the  elastic  release 
wave  crones  the  plastic  shock  near  the  flyer  free  surface  or 
target  free  surface,  a  discontinuity  of  0.01  g/cm  is  created 
in  iron.  When  a  rarefaction  wave  train  crosses  this  continuity, 
about  one  percent  of  it  is  reflected.  The  region  between  the 
spall  plane  and  the  free  surface  is  rich  in  such  discontinuities. 
Detailed  SWAP  calculations  show  that  the  characteristics  travers¬ 
ing  this  space  are  dispersed  and  attenuated.  The  consequence 
for  equation  (21)  is  that  spall  stress  will  be  underpredicted  by 
a  few  percent. 

The  second  source  of  errors  in  equation  (21)  is  overtaking 
of  the  SP-S  rarefaction  by  the  separation  shock.  Figure  9  shows 
the  effect  on  free  surface  motion;  the  solid  line  illustrates  the 
observed  free  surface  velocity  history,  and  the  dotted  line  shows 
the  history,  that  would  have  occurred  if  not  for  the  overtaking. 
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The  decrement  in  free  surface  velocity,  v'  is  lost,  and  cannot  be 

deduced  from  the  experimental  record.  Consequently,  equation  (21) 

will  underestimate  o  .  The  extent  to  which  spall  information  is 

s 

lost  depends  on  d,  and  whether  or  not  the  spall  shock  travels  at 
an  elastic  or  plastic  wavespeed.  All  right-travelling  release 
waves  arriving  at  the  spall  plane  for  a  time  A  .  before  spall 

3  U 

occurs  are  "lost".  For  T>^2d,  A  .  is  very  nearly 

S  t 


where  U  is  the  shock  velocity  from  SP,  and  UB  is  the  typical 

release  velocity.  It  has  also  been  assumed  that  SP  is  a  distance 

d  from  the  target  free  surface,  and  the  effects  of  changes  in 

particle  velocity  have  been  neglected.  For  steel,  A  ^0.1  us 

for  d  =  2mm  and  0.05  us  for  d  =  1  mm,  if  U  is  an  elastic  wave- 

s 

speed.  This  leads  to  a  substantial  value  of  v'  .  If  U  is  plas- 

s 

tic,  then  for  steel  At  v  0.016  us  for  d  =  2mm,  and  At.  v  0.008  us 

s  —  s  - 

for  d  =  1  mm.  SWAP  calculations  permit  an  estimate  of  the  magni¬ 
tude  of  these  corrections.  For  copper  and  other  targets  for 
which  the  elastic  strength  is  relatively  small,  the  corrections 
needed  for  this  effect  are  not  important. 

The  third  possibility  is  that  the  spall  plane  continues  to 

support  a  stress  after  fracture.  The  resulting  free  surface 

history  is  shown  in  Figure  10.  The  spall  signal  fails  to  rebound 

to  the  peak  free  surface  velocity,  missing  it  by  an  amount  Avf. 

Equation  (21)  is  still  valid  for  o  .  the  stress  at  which  the 

s 

fracture  initiates.  The  spall  plane  strength  is  given  by 

°1  -  ^  P0C8  4Vf  ,23) 

If  varies  with  time,  then  no  simple  interpretation  of  the 
velocity,  record  is  possible.  Most  spall  records,  at  least  near 
the  threshold,  look  like  Figure  10. 
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Figure  9.  Free  Surface  Velocity  History  Modified  by  Spall 
Shock  Overtake. 


s 


Figure  10.  Free  Surface  Velocity  History  Modified  by  Residual 
Stress  cn  Spall  Plane. 
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3.  THRESHOLD  CRITERIA 

The  data  obtained  in  this  program  were  analyzed  with  both 
time-dependent  and  time-independent  formulations.  The  Tuler- 
Butcher  relationship  was  used  for  the  time-dependent  approach 
(TULER  and  BUTCHER) .  The  mean  stress-plastic  strain  approach 
embedded  in  Figure  4  was  used  for  the  time-independent  approach. 

As  discussed  in  Section  1,  the  form  of  the  Tuler-Butcher 
equation  to  be  chiefly  employed  is 


K2  ‘  /<°  -  “o'"  dt 
K2  -  K2c  f°r  sPa11 


(24) 


In  equation  (24)  ,  oo  is  the  critical  stress  for  spall 

formation,  and  K  is  a  damage  parameter.  Spall  cannot  occur  if 

the  tensile  stress  does  not  exceed  o  .  The  tensile  loading  on 

the  incipient  spall  plane  is  approximately  a  pulse  of  amplitude 

-o  and  duration  2d/Cr,  (where  no  spall  occurs)  .  It  follows 
max  B 

that  the  threshold  conditions  for  formation  of  a  spall  plate  are 
approximately  as  shown  in  Figure  11a. 


Some  numerical  insight  into  the  consequences  of  equation  (24) 
can  be  obtained  by  substituting: 


0  *  1/2  po  uo  (CB  +  Suo> 
0  *  1/2  Co  CB  Uc 


(25) 


which  gives: 


2  K 

uo(l  *  1/2  Sucl  -  T  +  “<= 

U  o 


(26) 


Equation  (26)  is  graphed  in  Figure  lib  for  values  representative 
of  steels,  namely  p  =7.85  g/cm^ ,  C„  =  4.67  mm/ps,  S  =  1.33, u  = 
0.2  mm/ps,  and  K  =  337  kbar2^-  ps.  It  can  be  seen  from  Figure  lib 
that  the  range  of  flyer  plate  thickness  over  which  the  spall 
criterion  depends  on  thickness  is  rather  small,  d  <  1  mm.  The 
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Figure  11.  Spall  Threshold  Criterion  Predicted  by  Equation  (24) 
(a)  General;  (b)  Representation  of  Steel. 
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,  -1/2 
transition  from  d-dependence  to  uQ-dependence  depends  on  K  . 

It  is  necessary  to  use  impractically  large  values  of  d  to  observe 

spall  at  u  =  uc  directly.  Actual  points  on  the  curve  of  Figure 

11  can  be  used  to  calculate  oq  and  KXc.  For  such  a  threshold  shot, 

does  K.  depend  on  c  ,  ?  The  answer  can  be  obtained  by  using 
Ac  min  J  ^ 

SWAP  to  calculate  K,  for  shots  in  which  a  .  «  a  .  If  it  is 

A  min  o 

found  that  K ^  here  too,  then  the  answer  is  yes.  Otherwise, 

a  complete  spall  criterion  would  need  to  express  K  as  a  func- 

A  C 

tion  of  o  .  The  dependence  of  K,  on  a  can  have  two  origins, 
max  Ac  max  3 

Shock-induced  fractures  may  affect  K.  .  Or,  spall  may  be  asso- 
dated  with  a  critical  time  scale,  so  that  a  can  drop  far  below 
cq  for  a  while  before  spall  occurs.  This  latter  point  of  view 
is  more  consistent  with  recent  developments  in  fracture  mechanics. 

4.  VISAR  DATA  PRESENTATION 

Free  surface  velocities  were  measured  with  a  VISAR.  The 
analysis  of  VISAR  data  is  discussed  in  Appendix  C.  Each  VISAR 
record  yields  a  curve  showing  how  v  varies  with  time.  Uncer¬ 
tainties  in  VISAR  data  vary  from  shot  to  shot,  and  point  to  point. 
Time  of  arrival  data  are  almost  always  very  precise.  Likewise, 
relative  velocities  are  usually  rather  well  determined.  However, 
a  variety  of  effects  conspire  to  introduce  uncertainties  in 
absolute  velocity  measurements,  as  discussed  in  Appendix  C. 

In  order  to  express  uncertainties,  VISAR  results  are  often 
presented  in  tabular  form  in  this  report.  Figure  12  shows  a 
velocity  record  with  most  features  that  are  observed.  (Not 
every  shot  or  every  material  displays  all  of  these  features.) 

The  parameters  used  to  describe  the  VISAR  data  are  defined  in 
Figure  12. 
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Figure  12.  Representative  VISAR  Trace,  Showing  Parameter 
Definitions . 
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SECTION  III 

RESULTS  FOR  COPPER 


The  results  of  the  investigation  of  spall  in  copper  are 
organized  according  to  the  flow  chart  of  Figure  5. 

1  •  PREVIOUS  WORK 

Several  data  for  tensile  yielding  and  fracture  of  copper  can 
be  found  in  the  literature.  In  addition,  many  authors  have  tried 
to  use  ductile  fracture  in  copper  as  a  paradigm  for  fracture  of 
ductile  metals  in  general.  In  spite  of  these  previous  efforts,  a 
general  fracture  criterion  for  copper  for  use  in  HULL-type  finite 
difference  codes  does  not  seem  to  be  available. 

It  is  well  known  that  copper  is  a  rate  sensitive  material, 
and  there  is  evidence  that  the  phenomena  responsible  for  yielding 
and  fracture  change  at  strain  rates  below  about  100  s  1  (CHRISTMAN 
and  ISBELL,  1971;  BAUER  and  BLESS,  1978).  Therefore,  this  review 
will  be  mainly  limited  to  high  strain  rate  data. 

a .  Intermediate  Strain  Rate  Results 

CHRISTMAN  and  ISBELL  (1971)  performed  split  Hopkinson 
bar  tests  on  annealed  (RF35)  and  half-hard  (RF78)  OFHC  copper.  The 
results  of  uniaxial  stress  tests  are  shown  in  Figure  13.  Relative 
to  the  annealed  material,  the  half-hard  material  was  extremely  rate 
sensitive,  but  did  not  work  harden. 

Depending  on  hardness,  initial  yield  stress  varied 

from  0.3  to  3.3  kbar.  The  annealed  material  could  be  described  by 
work  hardening  exponents,  defined  in  the  equation: 


where  N  =  0.27.  For  half-hard  copper,  there  was  almost  no  work¬ 
hardening,  N  =  0.031.  The  exponential  work  hardening  relation  does 
not  fit  low  strain  rate  data  (PERRA,  1976). 
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TRUE  STRESS  (kbor) 


Figure  13. 


yielding  of  Half-Hard  and  Annealed 
One-Dimensional  Compressive  Stress 
and  ISBELL. 


OFHC  Copper  in 
from  CHRISTMAN 
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CHRISTMAN  and  ISBELL  also  investigated  the  extent  of 
the  Bauschinger  effect  in  copper  by  performing  compressive-tensile 
cyclic  loading  to  yield  at  low  strain  rates.  They  found  the  effect 
was  not  significant.  However,  MCMILLAN  et  al.,  1971,  found  it 
necessary  to  invoke  the  Bauschinger  effect  to  explain  shock  release 
data. 

SHULTZ  (1969)  determined  the  high  strain  rate  yield 

behavior  of  OFHC  copper  using  a  wire-impact  technique.  His  copper 

had  a  static  yield  strength  of  2.50  kbar,  comparable  to  the  half- 

3  -1 

hard  copper  tested  by  CHRISTMAN  and  ISBELL.  At  10  s  ,  this  material 
had  a  yield  strength  of  3.16  kbar,  which  is  consistent  with  CHRISTMAN 
and  ISBELL'S  data. 

Rapid  tensile  failure  in  a  uniaxial  stress  state  must 

be  accomplished  by  expanding  rings  or  cylinders.  These  have  been 

studied  by  BAUER  and  BLESS  (1979)  and  by  RAJENDRAN  and  FYFE  (1960). 

Both  sets  of  investigators  found  dramatic  increases  in  strain  to 

failure  due  to  suppression  of  local  instabilities  at  higher  strain 

rates.  This  increase  in  average  strain  to  failure  was  especially 

3  4-1 

large  in  the  interval  10  -10  s  .  According  to  BAUER  and  BLESS, 

ultimate  strain  (true  strain]  at  failure  actually  decreased  from 

.  •  3  -1 

a  static  value  of  1.6  to  value  of  1.0  at  e  ^  10  s  .  That  is, 

when  the  strain  in  the  failed  region  alone  is  considered,  it  becomes 

clear  that  the  material  is  less  ductile  at  high  rates,  as  shown  in 

Figure  14. 

b .  Spall  Data 

There  have  been  three  principal  experimental  investi¬ 
gations  of  spall  in  copper:  by  SRI,  summarized  in  SEAMAN  et  al.  1972; 
by  General  Motors,  summarized  by  CHRISTMAN  and  ISBELL,  1971;  and  by 
CEA,  reported  by  DAVID  et  al.  1975. 

SEAMAN  et  al.,  modelled  copper  with  the  SRI  NAG  theory. 
Parameters  in  the  theory  were  evaluated  by  impacts  below  the  spall 
separation  threshold.  Wedge  shaoed  flyers  were  sometimes  used,  in 
which  the  duration  of  the  tensile  pulse  varies  across  the  sample. 

The  material  studied  was  fully  annealed  OFHC  copper. 
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Eight  impacts  were  carried  out  at  velocities  up  to 
129  m/s,  which  was  below  the  separation  threshold.  The  void  con¬ 
centration  on  the  incipient  spall  planes  were  determined  by  counting 
on  sectioned  targets.  Several  parameters  in  the  NAG  theory  could 

not  be  directly  measured.  The  ii  itial  radius  of  voids  was  taken 

-4 

as  10  cm.  The  viscosity  of  copper  n  was  taken  as  750  p.  The 
constants  n  and  in  the  growth  aquation 

P-P0 

N  =  nQ  exp  — - -  (28) 

•  12  -3  -1 

.  _ _ _ _  _  o  o  ..  i  n1*1  - _ _  «*•  _ j  mu  — - 1. : 


were  nQ  =  2.8  x  10  cm  s  and  -  -2  kbar.  Then  the  nucleation 
threshold  turned  out  to  be  -5  kbar. 


In  relating  Pq  to  a  threshold  velocity,  the  yield 
behavior  of  the  copper  must  be  taken  into  account.  The  peak  tensile 
stress  in  planar  impacts  is  about  equal  to  the  peak  compressive 
stress.  For  a  given  strain,  the  peak  stress  in  one-dimensional 
strain  exceeds  the  mean  stress  by  2/3  Y.  If  pQ  =  5  kbar,  and  Y  'v 
3  kbar,  then  o  ^7  kbar,  which  means  that  spall  damage  will 

MclX 

start  to  accumulate  at  velocities  as  low  as  s,  25  m/s. 

The  growth  parameter,  PgQ,  defined  ins 


where  R  is  void  radius) ,  was  also  recommended  as  5  Kbar  for  copper. 

A  useful  equation  for  checking  the  NAG  formulation  is  the  prediction 
of  maximum  void  size  produced  by  a  square  tensile  pulse  of  duration 
At,  which  is: 


=  exp 


At 


In  the  SEAMAN  et  al.,  study,  the  NAG  parameters  were 
used  in  a  finite  difference  code  especially  modified  to  account  for 
the  influence  of  void  growth  on  mechanical  properties.  Naturally, 
this  formulation  had  to  account  in  an  average  way  for  the  plastic 
flow  and  plastic  instabilities  (see  PERRA,  1978)  that  occur  around 


voids.  The  agreement  with  predicted  void  concentrations  was  rather 
good,  but  void  content  was  underpredicted  at  high  damage  levels. 

This  indie, i  ^es  that  the  void  growth  assumption  employed  becomes 
invalid  as  voids  become  large  and  start  to  interact  with  each  other. 
SEAMAN  et  al.  did  not  present  an  explicit  discussion  of  uniqueness 
or  sensitivity  for  their  parameters,  although  they  did  discuss  sev¬ 
eral  test  calculations  which  also  yielded  reasonably  good  correla¬ 
tion  with  experiments.  In  these,  the  values  p^  =  6.5  kbar  and  g  = 
3570  p  gave  good  agreement. 

Doth  SEAMAN  et  al.,  and  MOSS,  (1977),  discuss  techniques 
for  predicting  NAG  parameters.  The  viscosity  parameter  that  deter¬ 
mines  void  growth  can  be  determined  independently.  WALTERS,  (1979), 
recently  reviewed  viscosity  data  and  reports  a  value  of  3000  p  for 
copper  at  93  kbar.  SEAMAN  et  al.  cite  literature  values  of  1150  to 
3600  p.  None  of  these  is  in  very  good  agreement  with  the  value  of 
750  p  needed  to  reproduce  the  observed  void  pattern  with  equation 
(29)  or  (30).  Similarly,  should  be  related  to  the  yield  stress 

and  fracture  toughness.  MOSS  derives: 


go 


K 


4R 


1C 


(31) 


However,  apparently  there  are  not  data  for  K^c  appropriate  for  shock- 
damaged  copper,  which  limits  the  usefulness  of  equation  (31). 

More  useful  is  the  identification  of  a  with  plastic 

go  r 

yielding  in  the  vicinity  of  voids.  CHRISTMAN  and  ISBELL  give  a 
value  for  the  HEL  of  8  kbar  (see  below).  SEAMAN  et  al.  identified 
void  growth  with  the  HEL  (p.57),  leading  to  a  prediction  of  pQ  = 

Orel  ”  2Y/3  which  agrees  with  the  value  of  p^  that  gave  gave  good 
results  in  calculations.  (It  is  nor  clear  in  their  discussion  if 
the  value  of  pn  was  in  fact  set  equal  to  oUT7T  -  2Y/3  and  not  varied)  . 

CHRISTMAN  and  ISBELL  determined  experimental  values  of 
spall  threshold  velocities  for  half-hard  OFHC  copper.  They  defined 
spall  threshold  as  the  point  at  which  cracking  occurred  over  50  per- 
■ent  of  the  spall  plane.  Their  data  are  give  in  Table  1.  Stress 
wave  profiles  were  measured  with  a  VISAR  or  other  techniques. 
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ULTIMATE  STRAIN 
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LOG  STRAIN  RATE  (S'1) 


Fiqu.ro  14.  Maximum  Average  Values  of  Strain  (Squares)  and  Local 
Values  of  Strain  (Circles)  as  a  Function  of  '  ‘•.rain 
Rate  for  Fully  Annealed  ETP  Copper.  From  BA>  f.r  and 
BLESS . 
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They  used  the  equation  (21)  to  calculate  the  spall  stress.  The 
result  depended  on  flyer  thickness,  being  o  =  -21  kbar  for  1-mm 

5 

flyers  and  -35  kbar  for  thinner  flyers. 

DAV1J  et  al.,  also  studied  spall  in  OFHC  copper  of 
0.01  mm  grain  diameter.  In  distinction  to  the  other  studies, 
they  used  impacP  velocities  in  the  range  2300  to  3500  m/s  to  pro¬ 
duce  shock  stresses  between  100  kbar  and  1  Mbar.  They  carried 
out  extensive  analysis  of  recovered  samples  and  observed  that 
cracks  occurred  in  parallel  zones  separated  typically  by  0.3  mm. 
(Some  of  the  results  c~  SEAMAN  et  al.  also  show  this.)  Measure¬ 
ments  ofAus,  accurate  only  to  about  40%,  are  given  in  Table  2. 

DAVID  et  al.,  propose  a  numerical  micromechanical  model 
quite  different  in  l pproach  from  that  of  SRI.  It  is  based  upon 
the  stress  needed  to  open  cracks  and  make  them  grow  in  essentially 
discontinuous  jumps.  The  critical  stress  for  growth,  o*,  is 
assumed  to  be  a  linear  function  of  pressure.  For  copper,  the 
model  parameters  are  <S  =  0 . 3  mm  (cracks  open  at  last  6/3  apart), 

0  =  2.55  us,  where  it  is  assumad  that  when  cracks  grow,  they  do 
so  at  a  rate: 


dc  c 
dt  "  0 


(32) 


c  being  crack  length,  and  the  stress  that  a  crack  can  support 


0  for  o  < 

2  2 

ql  for  c  > 

20a1 


(33) 


the  last  model  parameter  is  the  constant  k  defined  as: 

k  =  o*/P  (34) 

which  for  copper  is  equal  to  0.03.  Good  agreement  is  obtained  in 
finite  difference  calculations  of  free  surface  velocity  history 
using  this  rather  simple  model. 
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TABLE  2 

SPALL  RECOVERY  DATA,  FROM  DAVID  et  al. 


Flyer  Plate 
Velocity 
(m/s) 

Flyer 

Thickness 

(mm) 

Thickness 

(mm) 

Peak  Particle 
Velocity 
(m/s) 

AUg 

(m/s 

3500 

3 

20 

3200 

520 

25 

2800 

450 

30 

2530 

380 

35 

2235 

250 

40 

1190 

220 

2700 

3 

25 

2300 

250 

30 

2150 

230 

35 

1950 

200 

40 

1850 

120 

(80  mm  diameter  targets  and  flyers) 


Other  Properties  of  Co 


Table  3  contains  a  list  of  other  properties  of  OFHC 
copper  which  are  useful  in  finite  difference  modelling.  These 
values  are  taken  from  the  compilations  given  by  CHRISTMAN  and 
ISBELL  and  MCMILLAN  ec  al. 


2 .  EXPERIMENTAL  DETERMINATION  OF  SPALL  THRESHOLD 

The  copper  used  in  this  program  was  supplied  to  the  Univer¬ 
sity  by  AFATL.  It  was  in  the  form  of  a  bar,  33  mm  thick,  and 
it  was  labelled  as  cross-rolled  OFHC  plate.  Its  hardness  as  used 
in  these  impact  tests  was  RF48.  The  data  for  spall  threshold 
determination  is  shown  in  Table  4.  The  shots  used  copper  flyers 
and  targets,  except  shot  109,  in  which  an  aluminum  flyer  plate 
was  used  to  produce  lower  peak  pressures;  the  peak  pressure  in 
shot  109  is  equivalent  to  an  impact  with  a  copper  flyer  at  98  m/s. 
Peak  stress  values  were  calculated  from  standard  Hugoniot  data, 
namely,  U  =  3.9u  +  1.489  (in  units  of  mm/ys) . 

In  many  targets  extensive  void  production  was  observed, 
al chough  separation  at  the  spall  plane  did  not  occur.  Figure  15, 
for  example,  shows  a  section  through  the  target  in  shot  109,  and 
Figure  16  shows  a  view  of  the  spall  plane  viewed  on  a  surface 
parallel  to  the  spall  plane.  It  is  evident  that  the  mean  density 
on  this  plane  is  about  half  the  original  sample  density. 

Figure  17  shews  a  photomicrograph  from  shot  109  showing  damage  in 
the  trajectory  plane.  Figure  18  shows  a  macroscopic  view  of  damage 
caused  by  successively  higher  impact  velocities.  The  data  show 
that  for  the  1.5  mm  flyer  plates  against  4.0  mm  targets,  void 
growth  occurred  at  uQ  >_  244  m/s.  For  2  mm  flyer  plates,  spall 
separation  occurred  at  uq  >_  168  m/s.  This  is  significantly  less 
than  the  value  for  1.5  mm  plates  suggesting  rate  effects  are 
important.  Another  indication  of  results  were  obtained  in  shots 
189  and  322;  even  though  these  impacts  produced  the  same  peak 
stress,  the  results  differed.  For  2 . 5  mm  flyer  plates,  spall 
separation  occurred  at  286  m/s . 


TABLE  3 

PROPERTIES  OF  OFHC  COPPER 


Longitudinal  sound  speed  CT  =  4.757  mm/ys 

JL 

Shear  wave  speed  C  =  2.247  mm/ys 

.  s  3 

Density  p  =  8.917  g/cm 

Expansion  coefficient  B  =  51.5  x  10-^  +  .006  x  10-6T  °C-^ 

Heat  capacity  C  =  0.0915  +  0.00003T  cal/g 
P 

Melting  temperature  T  =  1080 °C 

m 

Boiling  temperature  =  2595°C 

Heat  of  fusion  Lf  =  49  cal/g 

Heat  of  vaporization  Ly  =  1146  ca.l  /g 

Bulk  modules  K  =  1418  kbar 
s 

9K  /9T  =  0.278T 
s 

9Kg/9P  =  5.53 

Shear  modules  y  =  450  kbar 
9y/9p  =  1.30 

Gruneisen  parameter  y  =  2.13 
9y/9T  =  0.00076  °C_1 

9y/9P  =  0.002  kbar-1 

.  2  2 
Low  pressure  Hugomot  =  1.5  +  3.44u  +  143u 

High  pressure  Hugoniot  (>1  Mbar)  U  =  3.964  +  1.463u 

Hugoniot  elastic  limit,  (anomalous  precursor  behavior,  best 
value) 


°HEL  =  8  kbar 


Release  wave  speed  C  >  5.3  mm/ys,  no  elastic/plastic 
separation  observed? 


TABLE  4 

COPPER  DATA  SHOTS 


Max 

Stress 

(Khar) 

Flyer 

Thickness 

(mm) 

— 

Target 

Thickness 

(mm) 

Shot 

No. 

Remarks 

80 

14 

1.50 

4.00 

121 

No  voids 

103 

19 

1.45 

3.90 

111 

Void  layer 

244 

45 

1.45 

4.00 

110 

Spall 

separation 

168 

30 

2.0 

4.00 

221 

Spall 

separation 

189 

35 

2.0 

4.00 

220 

Spall 

separation 

322 

35 

2.0  (Al) 

3.96 

109 

Void 

layer 

286 

53 

2.5 

3.94 

108 

Spall 

separation 

334 

62 

2.5 

4.00 

222 

Spall 

separation 

37 


-fcv.  ••  • 


1  mm 


Figure  16.  Photomicrograph  from  Shot  109,  Taken  in  Incipient 
Spall  Plane. 


0.1  mm 

Figure  17.  Photomicrograph  from  Shot  109,  Near  Separation  Plane. 


The  range  of  impact  parameters  over  which  void  formation 
occurs  is  comparatively  large  in  copper.  As  a  practical  manner, 
the  spall  criterion  used  by  CHRISTMAN  and  ISBELL  is  a  useful 
definition:  a  50  percent  reduction  in  cross  sectional  area  on 

the  incipient  spall  plane. 

3.  VISAE  DATA  FOP  COPPER 

Good  VISAR  data  were  obtained  on  two  copper  shots:  222  and 
220.  The  free  surface  records  were  distinctive  for  their  rather 
featureless  shape.  Figure  19  shows  the  oscilloscope  data  for 
shot  222.  The  "spall  signal"  is  very  small,  indicating  that  o 

5 

is  very  small.  The  VISAR  data  are  summarized  in  Table  5.  The 
parameters  in  Table  5  are  defined  in  Figure  12. 

4 .  THRESHOLD  CRITERIA 

The  threshold  data  from  this  work  and  from  CHRISTMAN  and 
ISBELL  are  plotted  in  Figure  20.  £'7AP  calculations  were  carried 

out  for  the  spall  threshold  data.  The  input  parameters  were 
based  on  the  conclusions  in  Section  III.l.  The  SWAP  model  for 
release  waves  was  a  simple  elastic  plastic  process  with  no  strain 
hardening  or  softening.  The  VISAR  data  did  not  reveal  sufficient 
structure  in  the  release  process  to  warrant  more  elaborate  models. 
The  parameter  was  calculated  from 

Kx  =  |(  |a-a0i)  Adt 

On  the  incipient  spall  plane,  taken  as  x  =  T  -  d.  The  nucleation 
stress  oQ  was  set  equal  to  -5  kbar,  based  on  the  results  discussed 
in  Section  3.1,  ar'd  \  =  2  as  discussed  in  the  Introduction.  The 
resulting  values  of  are  shown  in  Figure  21,  where  they  are 
plotted  against  flyer  plate  thickness  d. 

An  uncertainty  in  the  uc  data  of  CHRISTMAN  and  ISBELL  of 
10  m/s  should  be  applied  in  Figure  20  for  d  <  2  mm,  based  on  the 
scatter  in  their  results.  The  data  for  d  =  0.4  is  for  a  target 
in  which  the  peak  shock  stress  was  seriously  eroded  by  the  release 
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Figure  19.  VISAR  Record  for  Copper  Shot  222.  Fringe  Constant 
was  322  m/s. 


3  d  (mm) 


Figure  20.  Most  Tensile  Stress  vs  Flyer  Plate  Thickness  for 

Spall  in  Copper.  (Shading  ' ndicates  extent  of  spall 
Circles  are  present  work.  Flag  marks  shot  with 
aluminum  flyer  plate.) 


wave  from  the  flyer  rear  surface.  This  result  is  much  more  sen¬ 
sitive  to  the  treatment  of  the  release  process  in  SWAP  than  are 
the  other  data  points.  The  value  for  for  this  point  is  thus 
less  reliable. 

It  is  evident  from  Figure  21  that  the  value  of  is  a  much 

more  volatile  function  of  impact  velocity  than  is  o  .  The  data 

^  2  max 

can  be  described  by  50  <  K2c  <  110  kbar  -  us  with  two  exceptions. 

The  first  is  the  datum  at  d  =  0.4.  As  already  discussed,  this 

datum  is  considerably  less  precise  than  its  neighbors,  and  can 

probably  be  disregarded.  The  other  anomalous  point  is  from  shot 

109,  in  which  an  aluminum  flyer  was  launched  at  a  copper  target. 

In  the  SWAP  calculation,  the  description  of  aluminum  used  by 

BARKER  and  YOUNG  was  employed.  The  SWAP  calculations  of  K2c  are 

unlikely  to  be  in  serious  error.  Increasing  Oq  to  -10  kbar,  or 

decreasing  A  to  1  does  not  significantly  reduce  the  discrepancy 

between  this  shot  and  the  others  with  copper  flyer  plates.  The 

origin  of  the  discrepancy  has  not  been  resolved  at  this  time. 

Since  it  is  a  singular  occurrence,  it  is  recommended  that  a  value 

2 

of  K2c  =  50  to  110  kbar  -  ys  be  employed  as  a  spall  threshold 
predictor  for  copper. 

5.  OVERDRIVEN  SPALL 

The  spall  stress  computed  from  the  Av  data,  using  equation 

5 

(21)  turns  out  to  be: 


og  =  -3.2  t  0.8  kbar 

SWAP  calculations  also  require  values  of  ag  n,  -5  kbar  to  reproduce 
the  VISAR  data.  This  is  considerably  less  than  the  values  reported 
by  DAVID  et  al .  and  by  CHRISTMAN  and  ISBELL.  We  do  not  have  an 
explanation  for  why  the  values  of  o  observed  here  differ  from 
past  results.  We  note  that  the  ratio  of  d/T  was  only  half  that 
used  by  CHRISTMAN  and  ISBELL,  the  smaller  ratio  of  d/T  would 
probably  lead  to  a  large  value  of  v'  (as  discussed  in  Section  2). 
The  VISAR  data  also  showed  that  the  residual  stress  on  the  spall 
plate,  a^,  is  essentially  zero. 
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The  result  that  I  n  |  <  Jn^l  is  implausible.  One  way  to  re¬ 
concile  the  disparity  in  results  is  to  invoke  a  sudden  loss  in 
fringe  contrast  in  t  he  records  obtained  in  th  is  program.  CHRIST¬ 
MAN  and  I  KB F, LI.  used  a  specu lari  y-re  f  lee t  i nu  VISAR,  which  may  have 
been  less  sensitive  to  the  mechanism  responsible  for  contrast 
degradation.  on  the  other  hand,  loss  of  fringe  contrast  is  not: 
self-evident  in  the  data.  Perhaps  for  reasons  not  yet  determined, 
spall  signals  are  difficult  to  observe  in  copper.  It  is  recom¬ 
mended  that  this  problem  be  subjected  to  further  investigation. 

The  Av  data  of  CHRISTMAN  and  ISBELL  give  o  =  -20  t  1  kbar 
s  s 

for  throe  shots  with  1  mm  flyer  plates  and  o  =  -31  kbar  for  0.5 

mm  flyer  plates.  The  correction  to  these  values  due  to  the  Av ' 

effect:  discussed  in  Section  II  is  relatively  small,  since  copper 

has  relativo1y  low  yield  stress.  Calculations  of  K0  t  give  22  to 

53  kbar^  -  ijs  for  the  shot,,  with  d  =  1  mm,  and  24  kbar*"  -  ps  for 

the  shot  with  d  =  0.5  nun.  The  slight  discrepancy  between  these 

values  and  the  values  of  K„  derived  from  the  threshold  shots 

2c 

should  not  result  in  significant  errors  in  predictions  of  spall 

location  and  time  of  occurence.  For  the  present,  it  is  rccom- 

2 

mended  that  K2c  =  75  kbar  -  its  for  prediction  of  spall  of  OFHC 
copper . 


STRAIN  AT  FAILURE 


Between  threshold  conditions  and  overdriven  spall,  failure 

The  corres- 
022  . 


occurred  at  computed  values  of  -15  <  o  <  -31  kbar 


ponding  value  of  r.  ^ ,  from  SWAP  calculation,  are  0.01  <  c  < 

For  prediction  purposes,  a  mean  value  of  =  .015  can  be  used. 

The  value  of  o/Y  at  this  point  was  4.9.  The  uncertainty  in  this 

figure,  due  to  possible  errors  in  deviator ic  stress,  is  about 

±25  percent.  This  fracture  datum  should  be  compared  with  the 

datum  from  the  expanding  cylinder  experiments,  e  =  'vl.O  at 

P 

P/Y  =  0/3.  The  dependence  of  the  critical  value  of  on  P/Y  is 
very  significant.  However,  a  failure  criterion  by  itself 
misses  the  essential  time-dependence  of  spall-type  fracture  in 


copper , 
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SECTION  IV 

RESULTS  FOR  SAE  1020  STEEL 


1.  PREVIOUS  WORK 

The  mechanical  properties  of  mild  steel  have  been  discussed 
by  many  authors.  Recent  data,  as  well  as  a  pertinent  review,  can 
be  found  in  BLESS  and  BARBER  (1979).  In  their  study,  the  1020 
steel  hardness  was  Rockwell  B  81*2.  Tensile  tests  were  performed 
with  conventional  testimj  machines  and  with  a  split  Hopkinson  bar 
device  operated  by  Dr.  T.  Nicholas  at  the  Materials  Laboratory, 

Air  Force  Wright  Aeronautical  Laboratories.  The  material  was 
found  to  be  significantly  strain-rate  dependent.  The  yield 

strength  associated  with  the  recovery  from  the  overshoot  increased 

-4  1 

from  3.29  kbar  at  10  /s  to  4.1  kbar  at  1/s,  to  6.6  kbar  at  10  /s. 

There  was  also  no  downward  turning  of  the  stress  strain  curve  at 

the  highest  strain  rate;  this  is  apparently  associated  with  the 

3 

relative  decrease  in  necking  with  strain  rate.  At  10  /s  the 
sample  could  not  be  driven  to  rupture  because  the  load  duration 
was  limited  to  300  p s .  The  true  strain  at  failure  was  estimated 
from  the  cross  sectional  area  of  the  failed  region  to  be  0.9. 

The  SRI  uucleation  and  growth  model  for  fracture  has  been 
applied  to  Armco  iron.  It  is  a  priori  plausible  that  these  param¬ 
eters  will  apply  to  mild  steel  also.  As  given  by  SEAMAN  and 

SiiOCKEY  (  1972)  ,  the  nucleation  threshold  stress  is  3  kbar  (p  in 

o 

equation  28).  The  sensitivity  parameter,  p^  is  4.56  kbar.  The 
growth  threshold,  p  in  equation  (29),  is  -0.2  kbar. 

Matuska  has  carried  out  Hull  code  calculations  to  try  to 

reproduce  the  hemispherical  impact  data  of  BLESS  and  BARBER. 

He  employed  a  work  hardening  model  designed  to  fit  the  observed 

velocity  decay  -  ¥  =  4.6  kbar,  increasing  linearly  to  6.3  kbar 

at  '  -  0.3,  and  thereafter  constant.  He  used  an  ultimate  failure 

criteria  discussed  in  Section  1.2.  The  ultimate  failure  surface 

was  fit  to  the  Hopkinson  bar  results.  It  passed  through  the 

points  P/Y  =  0.33  and  t  =  0.9,  and  P/Y  -  1.5  i  =  0.3.  The  HULL 

P  P 
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Figure  22 


Comparison  of  Measured  and  C^lc^^®dBfpg|SSB^pBEp 
in  Impacted  1020  Steel  Targets  (from  BLESS,  BARBER, 

and  MATUSKA) . 
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calculations  were  able  to  reproduce  the  observed  velocity  decay  and 
failure  patterns.  This  agreement  was  only  obtained  after  inter¬ 
actions  on  yield  and  failure  criteria.  Figure  23  shows  the  calcu¬ 
lated  failure  zones  for  the  experiment  pictured  in  Figure  1.  This 
result  is  described  by  BLESS,  BARBER,  and  MATUSKA.  These  fracture 
data  for  1020  steel  are  graphed  in  Figure  4a. 

The  variation  of  release  wave  speed  and  bulk  sound  speed 
with  pressure  is  essentially  linear  for  a-ferrite.  This  is  dis¬ 
cussed  in  Appendix  B. 

2.  EXPERIMENTAL  DETERMINATION  OF  SPALL  THRESHOLD 

The  material  used  in  the  study  was  commercial  SAE  1020  steel. 
Representative  microstructures  are  shown  in  Figure  24.  The  hard¬ 
ness  varied  from  sample  to  sample.  The  range  for  almost  all 
targets  was  RB  62  to  75. 

The  data  shots  for  spall  threshold  for  1020  steel  are  sum¬ 
marized  in  Table  6.  The  compressive  stress  is  calculated  from 
the  a-phase  iron  Hugoniot  (VAN  THIEL  et  al.,  1968). 

U  =  4.62  +  1.74  u 

Tensile  stress  maxima  were  computed  with  the  SWAP  code,  as 
described  below.  There  was  no  discernible  variation  in  the  position 
of  the  spall  plane  for  2  mm  flyers;  the  values  were  all  nearly 
1.9  mm.  In  the  shots  with  1  mm  flyers,  spall  thickness  was  only 
0.65  mm  to  0,79  mm.  The  experimental  spall  thresholds  may  be 
summarized  as  follows:  for  1  mm  flyers  and  5  mm  targets, 

251  <  u  <  346  m/s.  For  2  mm  flyers  and  5  mm  targets,  a  v  :'.J9  m/s. 

U  V 

Measurements  were  made  of  m.icrohardness  in  shot  58  across  a 
cross  section.  In  tht.  center  of  the  shocked  region  the  hardness 
was  VH  148.  No  softening  on  the  incipient  spall  plane  was  detected. 

The  evidence  from  the  threshold  experiments  indicates  that 
under  the  circumstances  considered  here,  tensile  rupture  is  so 
abrupt  that  it  would  probably  not  be  useful  to  consider  a  partial- 
damage  model.  In  only  two  shots  were  voids  nucleated  that  failed 
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Figure  23.  Agreement  to  Observation  of  Figure  1  Obtained  With 
Hull  Code  Calculation  (From  BLESS,  BARBER  and 
MATUSKA) . 
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SUMMARY  OF  SPALL  OBSERVATIONS  FOR  1020  STEEL 
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to  coalesce  into  a  rupture  plane.  Shots  differing  in  ^max  by  only 
4  kbar  showed  complete  spall  or  no  spall.  Figure  25  illustrates 
the  appearance  of  voids  in  a  sub-threshold  shot. 

3 .  FREE  SURFACE  VELOCITY 

VISAR  data  were  obtained  for  both  loading  and  unloading  of 
this  material.  A  summary  of  the  shot  configurations  for  which 
VISAR  data  were  obtained  is  given  in  Table  7.  The  free  surface 
trace  parameters  are  given  in  Table  8.  The  elastic  precursor  was 
resolved  in  three  shots  as  shown  in  Figure  26.  The  stress  levels 
corresponding  to  these  shots  are  approximately  9,  27,  and  42  kbar, 
all  for  4  mm  thick  targets.  Only  the  low  value  is  consistent  with 
HEL  data  for  a-iron:  BARKER  and  HOLLENBACK  found  10.5  kbar  for 
6.4  mm  thick  samples;  MCQUEEN  et  al.  give  14  to  15  kbar  for  "hard' 
iron,  and  9  to  11  kbar  for  "soft"  iron.  The  observed  higher 
amplitude  arrivals  may  be  reflections  of  the  elastic  wave  from 
the  advancing  plastic  wave  front;  the  data  for  both  shots  14  and 
23  are  consistent  with  indistinct  arrivals  at  ^14  kbar.  We  con¬ 
clude  that  the  initial  elastic  wave  was  missed  in  those  shots. 

Figure  27  shows  sample  curves  for  release  wave  arrivals  and 
spall  signals.  These  shots  were  very  helpful  in  developing  a 
deviatoric  stress  model  and  spall  criteria  for  use  in  SWAP  cal¬ 
culations,  as  explained  below. 

4 .  SWAP  MODEL 

The  SWAP  models  were  tried  using  various  descriptions  of 
unloading,  given  in  Table  9.  All  models  employed  identical  loading 
description,  appropriate  for  amax  <  130  kbar,  namely: 

3 

pQ  =  7.846  g/cm 

CQ  =  3.61  mm/ys  (VAN  THIEL,  3  968) 

S  =  1.75  (VAN  THIEL,  1968) 
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Figure  26. 
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Observations  of  Elastic  Precursor  Signals. 
14  (Closed  Circles)  63  (Open  Circles)  ,  and 
(Squares) 
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Figure  27 


VISAR  DATA  from  Shots  24  and  63 
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VISAR  DATA  AVAILABLE  FOR  1020  TARGETS 


Shot 

d 

(nun) 

T 

(mm) 

Qual: 

223 

1.0 

4.0 

Good 

63 

1.0 

4.0 

Good 

14 

2.0 

4.0 

Good 

24 

2.0 

5.0 

Good 

13 

2.0 

4.0 

Good 

273 

2.0 

Cone 

Fair 

274 

2.0 

Cone 

Good 

275 

2.0 

Cone 

Good 
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TABLE  8 


FREE  SURFACE  VELOCITY  PARAMETERS  FOR  1020  STEEL 
(Units  =  ys  and  m/s) 


SHOT 

223 

63 

14 

24 

13 

.15  i  .03 

.18  ±  .02 

.16  ±  .02 

.10  +  .02 

.12  ±  .05 

t 

.22  +  .01 

.23  ±  .01 

.67  ±  .02 

.16  ±  .02 

.28  ±  .02 

m 

t  , 
rl 

.40  ±  .01 

.24  ±  .01 

.78  ±  .02 

.61  ±  .01 

.65  ±  .02 

fcr2 

.43  ±  .01 

.31  ±  .01 

NR 

NR 

NR 

tr3 

NR 

.46  ±  .01 

NR 

NR 

NR 

fcs 

.53  +  .01 

.51  ±  .01 

.  88? 

.90  ±  .02 

.85  ±  .02 

tB 

.57  +  .01 

.56  ±  .01 

NR 

.97  ±  .02 

.92  +  .02 

V 

576  ±  15 

574  ±  10 

950  ±  15 

600  i  15 

550  +  10 

p 

V  . 

940  ±  10 

428  •  10 

685  +  15 

490  ±  10 

430  ±  5 

min 

VE 

115  ±  10 

39  ±  5 

180  ±  10 

NR 

Av;. 

140  t  10 

146  ±  10 

265? 

100  ±  10 

120  ±  10 

Avf 

20  ±  5 

92  +  10 

NR 

6  5  +  10 

90  1  15 

Vr2 

25  +  5 

16  ±  5 

NR 

NR 

NR 

vr3 

NR 

81  +  10 

NR 

NR 

NR 
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TABLE  9 


VALUES  OF  YIELD  PARAMETERS  FOR  VARIOUS  SWAP  MODELS 
Model  C  gave  the  best  agreement  with  VISAR  data. 


Model 

A 

B 

C 

Y  eq . 

1 

2 

3 

Y 

o 

7.124 

6.50 

6.30 

yl 

0.3125 

3.20 

4.60 

*2 

7.36 

60.68 

0.3 

W 

y  eq.  1:  Y  =  Y  (l+y.,c)(l — rr-)  (W  is  internal  energy) 

°  1  x2 

y  eq.  2:  Y  =  YQ(lfaP)  (P  is  hydrostatic  stress) 

a  =  (y^/yQ-l) /Y2  y^  is  max  Y  at  P  =  y2 

(Cp  is  plastic  strain) 
y^  is  max  Y  at  ^  =  y2 


y  eq .  3:  Y  =  Y  (1+ut.  ) 
op 


Y  =  1.25  (BARKER  and  YOUNG,  1974) 
dy/de  =  -10  (BARKER  and  YOUNG,  1974) 

3y/9P  =  0 

Model  A  has  a  constant  yield  stress  that  has  been  adjusted  to 
describe  our  observed  elastic  processor.  Model  B  includes  the 
yield  description  used  by  BARKER  and  HOLENBACK  1972,  for  Armco 
iron.  Model  C  (strain  hardening)  is  the  one  used  by  Matuska  to 
duplicate  the  results  of  BLESS  and  BARBER.  The  differences  in 
assumed  yield  stress  cause  the  variations  in  the  amplitudes  of 
the  first  two  release  arrivals  which  are  elastic  reflections 
from  the  elastic  precursor  and  main  shock  at  the  back  of  the  flyer 
plate.  Figure  28,  for  example,  illustrates  the  release  wave  fronts 
generated  by  these  models,  crossing  the  impact  plane.  Model  C 
clearly  gave  the  best  agreement  with  data.  Figure  29  shows  the 
SWAP  calculations  for  shot  24.  The  agreement  with  the  VISAR 
data  is  good.  Agreement  in  shot  63  was  equally  good.  Shot  223 
(Figure  30),  however,  showed  the  elastic  wave  arriving  a  little 
faster  than  in  the  SWAP  model.  It  would  appear  that  only  a 
varying  shear  modulus  could  reproduce  this  result.  In  spite  of 
this  defect,  the  predicted  time  of  arrival  and  amplitude  of  the 
spall  signal  are  in  excellent  agreement  with  the  data.  Overall, 
the  SWAP  model  was  judged  adequate  for  computing  stress  histories 
on  spall  planes. 

5.  SPALL  THRESHOLD  CRITERIA 

Spall  threshold  data  are  graphically  summarized  in  Figure  31. 
The  SRI  data  for  Armco  iron  are  also  included.  It  can  be  seen 
that  the  1020  steel  tested  here  is  clearly  more  resistant  to 
fracture  than  the  iron. 

SWAP  runs  were  carried  out  to  calculate  Kg  on  the  incipient 
spall  planes  for  1  mm  and  2  mm  flyer  plate  impacts.  The  parameter 

was  varied  until  the  same  value  of  K~  was  obtained  for  both 
o  2c 

conditions.  The  optimum  value  was  ao  =  -29  ±  1  kbar.  For  aQ  = 

-29  kbar,  spall  occurred  when  Kg  =  34  ±  1  kbar^  -ys. 


60 


*  0.6  065  0.7  Q75  Q8 

TIME  AFTER  IMPACT  (^s> 


Figure  28.  Release  Waves  from  Flyer  Free  Surface  Crossing 
the  Impact  Plane  in  Calculation  for  Shot  24 


TIME  AFTER  IMPACT  {fi  s) 

Figure  29.  SWAP  Result  for  Shot  24. 
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Figure  31. 


Data  for  Spall  Threshold  in  1020  Steel. 
(Squares  are  data  from  SEAMAN  et  al.; 
relative  shading  indicates  extent  of  spall.) 


63 


-  v,  « 


OVERDRIVEN  SPALL 


6  . 

The  spall  stress  calculated  from  the  VISAR  data  using  equation 

(21)  varies  with  flyer  thickness.  For  d  =  1  mm,  a  =  -34  ±  1  kbar. 

s 

For  d  =  2  mm,  o =  -25  ±  3  kbar.  This  is  consistent  with  the  pre¬ 
sent  damage  model.  The  release  fan  from  thinner  flyers  is  tighter, 
and  thus  the  value  of  oQ  may  be  more  exceeded. 

Attempts  to  use  SWAP  to  reproduce  observed  spall  signals  were 
not  wholly  satisfactory.  The  magnitude  of  Avg  was  relatively  insen¬ 
sitive  to  a  .  Values  -40  to  -50  kbar  provided  the  best  agreement, 
s 

but  the  uncertainty  is  about  ±5  kbar.  The  problem  seemed  to  be 
that  even  finely  zoned  problems  digitized  release  waves  in  'v-  10 
kbar  steps.  The  SWAP  calculations  showed  the  *~pall  shock  travel¬ 
ling  at  the  longitudinal  elastic  sound  velocity,  and  that  value 

was  used  to  calculate  the  values  of  a  given  in  the  preceding  para- 

s 

graph.  The  SWAP  calculations  showed  that  the  errors  in  using 
equation  (21)  may  amount  to  as  much  as  a  factor  two  in  spall  stress. 
It  is,  therefore,  concluded  that  the  most  accurate  measures  of  a 
are  the  relatively  imprecise  values  required  to  bring  SWAP  and 
VISAR  data  into  agreement.  In  other  words,  cg  =  -45  ±  5  kbar  for 
both  1  mm  and  2  mm  thick  flyers,  when  u  »  uc> 

The  SWAP  code  was  used  to  estimate  K2c  in  cases  of  rapid  spall. 

It  turned  out  that  the  values  were  very  small,  between  5  and  10 
2 

kbar  -  us.  Uncertainty  in  K2c  is  probably  of  little  practical 
importance  because  exact  values  of  a  and  K-  have  very  little 
effect  on  the  location  or  time  of  spall  fracture.  In  most  impact 
situations,  when  K2  is  exceeded,  it  is  exceeded  very  quickly,  and 
it  is  first  exceeded  in  the  region  where  3o/3t  has  the  most  nega¬ 
tive  value.  A  more  complete  model  is  obtained  by  considering  K2c 

a  function  of  a  .  ,  as  shown  in  Figure  32.  According  to  this 
mm 

model,  the  threshold,  criterion  involving  K2c  and  is  only  impor¬ 
tant  near  the  threshold,  where  a  .  does  not  exceed  a  by  more 

mm  o 

than  about  10  kbar.  For  more  extreme  value  of  a  .  ,  spall  takes 

mm 

place  more  rapidly  than  would  be  predicted  by  a  constant  value  of 

K-  .  This  conclusion  is  relatively  insensitive  to  the  value  of  a  . 
2c  1  s 
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7. 


TWO-DIMENSIONAL  TARGETS 


A  limited  number  of  experiments  were  done  with  two-dimensional 
impacts  in  order  to  more  effectively  assess  the  effects  of  shear 
strain  on  fracture.  The  targets  were  cones  25  mm  in  height  with 
a  90°  included  angle.  They  were  machined  from  1020  steel.  The 
flyer  plates  were  2  mm  thick,  38  mm  diameter,  1020  steel  discs. 

The  shot  matrix  is  presented  in  Table  10.  Predictions  of  the  re¬ 
sults  of  the  impact  were  made  by  staff  at  the  Air  Force  Weapons 
Laboratory  (AFWAL)  using  the  HULL  code.  Figure  33  illustrates  a 
target  before  impact.  A  single  crush  pin  was  emplaced  to  provide 
a  trigger  for  instrumentation.  Table  10  lists  the  shots. 

a .  Fracture  of  Conical  Targets 

Figure  34  through  36  show  cones  recovered  from  several 
shots.  It  can  be  seen  that  fracture  developed  on  four  separate  sys¬ 
tems.  A  spall- type  fracture  developed  parallel  to  the  cone  surfaces. 
This  fracture  was  apparently  the  first  one  to  develop,  and  it  is 
first  manifest  at  about  midheight. 

The  second  fracture  system  consisted  of  radial  cracks 
which  initiated  on  the  cone  surface.  These  cracks  seemed  to  form 
at  separation  angles  of  35  to  40°.  They  were  influenced  by  the 
presence  of  the  crush  pin  hole.  The  jagged  appearance  of  the  cracks 
suggests  that  they  were  formed  very  rapidly  over  much  of  the  surface 
of  the  cone.  The  earliest  cracks  coalesced,  and  the  later  cracks 
were  arrested. 

The  third  crack  system  was  a  roughly  cylindrical 
failure  abouc  12  mm  in  diameter  symmetric  to  the  cone  axis.  This 
failure  resembles  base  failures  in  impacting  cylinders  observed  by 
PAPIRNO  et  al.  Lastly,  there  was  a  radial  crack  system  expending 
outward  from  the  cylindrical  fracture  surface.  Its  trace  was  clearly 
visible  on  the  impact  surface.  These  cracks  were  spaced  by  about 
90°  . 

b .  Surface  Velocity 

Surface  velocity  was  measured  by  Hopkinson  flyoff 
plates  and  by  the  VISAR.  Hopkinson  flyoff  plates  are  small  discs 


Figure  33.  Preimpact  Photograph  of  Conical  Target 
from  Shot  276  . 


Figure  34.  Cross  Section  of  Recovered  Cone  in  Shot  273, 
u  =  383  m/s. 
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Ficiure  37. 


Flyoff  Plates  in  Motion  from  Shot  277. 


Figure  38. 


Superimposition  of  Preinpact  and  Postimpact 
Frames  from  Shot  278.  (Separated  by  42  jjs) 


that:  are  placed  on  the  target  surf  a  -e ;  they  are  launched  at  the 
peak  free  surface  velocity.  In  shot  277,  the  discs  were  0.05  mm 
steel  shims,  1.7  mm  in  diameter  (see  Figure  33).  They  were  con¬ 
toured  to  the  target  surface  and  held  in  place  by  minute  amounts 
of  silicon  greaS'  .  In  shot  278,  the  disc  size  was  decreased  to 
1.2  mm,  and  they  were  held  in  place  by  magnetization.  The  flyoff 
plates  were  placed  on  the  two  sides  of  the  target  silhouette.  The 
motion  of  the  flyoff  plates  was  photographed  by  a  Beckman  and 
Whitley  Model  300  framing  camera  at  a  rate  of  about  1.2  x  10 
frames/s.  Sample  frames  from  shot  277  are  shown  in  Figure  37  and 
38. 

Results  of  limited  precision  were  obtained  from  the 
two  flyoff  plate  shots.  In  shot  277  flyoff  plates  were  relatively 
large  and  they  only  moved  about  100  mm  during  the  camera  recording 
time.  The  normal  velocity  near  the  cone  apex  was  260  +  20  m/s.  At 
midhcight  it  was  310  '  30  m/s,  and  near  the  base  it  was  290  *  30  m/s 
The  flyoff  plates  in  shot  278  turned  out  to  be  too  small  to  reliably 
observe,  and  no  velocity  data  were  obtained.  The  photographic  data 
does  clearly  show  the  bulging  of  the  target.  The  target  reaches  its 
final  shape  by  40  us  after  impact.  Reliability  of  flyoff  plate  data 
should  be  checked  by  using  plates  of  varying  size  in  order  to  show 
that  the  size  of  the  geometry  disc  does  not  affect  launch  velocities 
Since  this  test  was  no<  done,  the  quoted  results  should  be  regarded 
as  upper  bounds. 

Three  shots  were  carried  out  using  the  VISAR  to 
measure  free  surface  velocity  history  at  a  point  7  ram  from  the 
cone  base.  In  one,  only  a  digital  oscilloscope  recorded  the  data 
(due  to  triggering  error,  and  the  digitizing  rate  of  50  ns  was  too 
low  to  reliably  follow  the  signal.  Oscilloscope  data  from  the  other 
two  shots  are  shown  in  Figures  39  and  40. 

The  normal  velocity  history  is  shown  in  Figure  41. 

It  has  been  assumed  there  that  no  fringes  were  lost.  The  data  are 
ambiguous  in  that  the  velocity  at  all  points  beyond  the  break  in 
the  curve  could  exceed  that  shown  in  the  figure  by  322  m/s. 
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Figure  39.  VISAR  Record  of  Normal  Surface  Velocity  from 
Shot  274  . 


Figure  40.  VISAR  Record  of  Shear  Surface  Velocity  from 
Shot  275  . 
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FREE  SURFACE  VELOCITY 
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The  uncertainty  in  velocity  can  be  resolved  with 
the  aid  of  the  data  from  shot  275.  In  this  show  the  VISAR  was 
incident  at  an  angle  of  45°,  and  viewed  at  an  angle  of  45°. 

The  velocity  measured  by  the  VISAR  is  equal  to  2vsin0^,  where 
v  is  the  free  surface  normal  velocity  and  0^  is  the  angle  of 
incidence.  (Initially  an  error  was  made  in  this  derivation 
and  it  was  thought  that  the  measured  velocity  was  sensitive 
to  the  in-plane  motion  of  the  free  surface;  measurement  of  the 
in-plane  motion  was  the  motivation  for  the  shot) .  Due  to  the 
spread  of  viewing  angles  (about  ±15°),  there  is  probably  loss 
of  fringe  contrast  in  shot  275.  The  peak  measured  velocity 
in  shot  275,  uncorrected  for  change  in  fringe  contrast,  is 
59  m/s,  giving  V  =  83  m/s.  The  actual  velocity  could  have  been 
higher  by  about  100  m/s,  due  to  loss  of  fringe  contrast.  If 
the  peak  normal  surface  velocity  was  actually  522  m/s  (corres¬ 
ponding  to  loss  of  a  f .  inge  in  shot  275),  then  the  measured 
peak  velocity  in  shot  275  would  have  been  369  m/s,  which  means 
more  than  one  fringe  would  have  occurred  in  the  first  half 
microsecond.  This  interpretation  seems  clearly  incompatible 
with  the  oscillograph.  Thus,  the  data  from  shot  276  support 
the  interpretation  given  in  Figure  41. 

c •  Shock  Arrival  Velocity 

The  shock  arrival  velocity  along  the  cone  free 
surface  was  measured  in  shot  276.  A  rear-surfaced  mirror  was 
cemented  to  the  surface.  The  rupture  of  the  mirror  was 
observed  using  a  pulsed  light  source  and  a  70  mm  streak  camera. 
The  velocity  14  mm  above  the  base  was  5.0  ±  .05  mm/ys.  20  mm 
above  the  base  the  velocity  was  4.61  ±  .14  mm/ys.  These  are 
consistent  with  the  plastic  wave  speed  in  steel. 

d .  Code  Simulation 

A  HULL  code  simulation  of  the  cone  impacts  was 
carried  out  by  personnel  from  AFWL.  The  impact  velocity  was 
500  m/s.  No  fracture  criterion  was  used.  Instead,  the  density 
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Figure  42.  Density  Contours  3.548  vis  After  Impact,  from  a 
HULL  Calculation,  Showing  Regions  Where  p  <  p 
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was  allowed  to  drop  to  less  than  p  .  The  result  is  shown  in 
Figure  42  for  the  time  at  the  end  of  the  calculation. 

The  failure  parallel  to  the  free  surface  has  been 
predicted.  However,  the  fracture  takes  place  1.5  mm  from  the 
surface,  whereas,  actually  the  distance  was  3  mm.  The  radial 
fracture  patterns  and  interior  conical  fractures  have  not 
developed.  At  this  time  there  are  no  negative  stresses  ex¬ 
ceeding  -15  kb,  so  further  fractures  near  the  cone  base  are 
unlikely.  However,  the  compressive  wave  has  not  yet  reached 
the  apex,  so  additional  fracture  may  occur  in  the  upper  part 
of  the  cone.  In  conclusion,  it  appears  that  the  fracture 
model  in  the  code  is  not  adequate  to  predict  at  least  some 
of  the  major  features  of  the  target  response. 

The  peak  free  surface  velocity  predicted  by  the 
calculation  is  200  m/s  at  15  mm  height,  and  460  m/s  at  7  mm 
height.  These  differ  by  about  50%  from  the  measured  values. 
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SECTION  V 

RESULTS  FOR  4340  STEEL 

1 .  LITERATURE  REVIEW 

Most  of  the  physical  properties  of  AISI  4340  steel  do  not 
differ  significantly  from  SAE  1020  steel.  The  low  pressure 
shock  properties  have  recently  been  reported  by  FRANZ  and 
ROBITAILLE.  The  a-phase  Hugoniot  is:  U  =  4,67  +  1.44u  (for 
units  of  mm/us).  According  to  these  authors,  the  a  -*■  e  phase 
transition  takes  place  at  132  kbar.  The  longitudinal  elastic 
wave  speed  is  5.86  mm/us,  and  the  bulk  wave  speed  is  4.67  mm/us. 
The  e  phase  Hugoniot  of  4340  is  the  same  as  pure  iron  within 
experimental  error. 

The  Hugoniot  Elastic  Limit  (HEL)  of  4340  steel  varies 
widely.  GUST  reports  apparently  random  variations  from  10  to 
28  kbar  in  RC58  material.  Apparently,  variation  of  10%  in  c>H  L 
may  be  due  to  effects  of  sample  thickness.  (The  HEL  decreases 
with  propagation  distance  due  to  dissipation) .  GUST  summarized 
HEL  data.  A  p] ot  cf  his  summary,  together  with  data  of  FRANZ 
and  ROEITAILLE  is  shown  in  Figure  43.  For  the  material  used 
in  this  study,  oUT;,T  =  20  kbar  is  appropriate. 

Hhiij 

UAPIRNO  et  al.  found  that  for  RC  52  material,  the  dynamic 
yield  point  was  the  same  as  the  static  value  at  50  percent  strain, 
namely  25  kbar.  BLESS  and  BARBER  report  dynamic  stress  strain 
curves  of  Nicholas  for  annealed  material.  Neglecting  a  short 
"overshoot",  characteristic  of  high  strain  rate  tests,  the 
yield  strength  increased  from  8.2  kbar  to  11.3  kbar  as  l  increased 
from  10‘V1  to  103s~^.  The  ultimate  strain  was  0.16.  BUTCHER 
and  CANON  also  studied  the  effect  of  strain  rate  on  yield.  Soft 
material,  RC  15  to  RC  32  was  quite  rate  dependent;  however, 

RC  54  material  was  not. 

There  do  not  appear  to  have  been  previous  studies  of  spall 
in  4340  similar  to  the  present  study.  Cylinder  impacts  of  4340 
were  carried  out  by  MESCALL  and  PAPIRNC,  Their  tests  were  with 
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RC15  material.  They  used  a  simple  wave  propagation  code  and 

achieved  satisfactory  agreement  between  post-shot  observation 

of  tensile  failure  and  calculations.  The  spall  model  of  a 

Tuler-Butcher  form  was  used  with  X  =  2,  o  =7  kbar,  and  K„  = 

2  o  2c 

337  kbar  y  s .  These  results  were  probably  not  very  sensitive 
to  the  spall  threshold  criterion,  since  peak  stresses  were 
well  above  uq.  NOVIKOV,  et  al.  (1966),  studied  a  Russian  steel 
which  is  probably  similar  to  4340.  In  the  annealed  version 
(Y  =  4.4  kbar),  the  critical  rupture  stress  was  28.4  kbar. 

Rupture  stress  was  calculated  from  free  surface  history  measure¬ 
ments  by  equation  (21) .  Thus,  the  rupture  stress  is  similar 
to  PSPALL  as  used  in  SWAP  calculations. 

Static  studies  of  tensile  fracture  in  4340  have  been 
carried  out  by  several  authors.  COX  and  LOW  reported  valuable 
observations  and  summarized  previous  findings.  Tensile  fail¬ 
ure  initiates  with  spherical  voids  which  grow  at  inclusion  sites. 
The  larger  inclusions,  principally  MnS,  nucleate  voids  first. 

The  higher  the  purity,  the  higher  the  stress  required  to  achieve 
a  given  void  content.  Voids  coalesce  by  propagation  of  hairline 
cracks  which  involve  fracture  initiation  at  tiny  cementite 
inclusions.  The  first  hairline  cracks  grow  parallel  to  the 
maximum  shear  stress.  The  amount  of  triaxiality  in  the  stress 
significantly  effected  void  growth  and  coalescence. 

COX  and  LOW  measured  the  fraction  of  void  area  on  the 
cross  sections.  This  was  found  to  vary  with  stress  state  and 
material  purity.  It  did  not  correlate  with  fracture  strain. 

A  typical  curve  is  shown  in  Figure  44. 

2.  DATA  FOR  SPALL  THRESHOLD 

The  4340  material  tested  in  this  program  was  the  same  as 
used  in  BLESS  and  BARBER.  It  was  vacuum  remelted,  and  was 
RC30  ±  2. 
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The  data  for  spall  threshold  for  4340  steel  are  shown 
in  Table  11.  For  1  mm  flyers,  there  were  some  targets  for 
which  partial  spall  damage  was  observed.  The  range  of  stress 
over  which  this  occurs  was  at  least  6  kbar.  In  shot  56,  although 
some  voids  were  present  in  the  polished  section,  the  VISAR 
records  did  not  indicate  a  spall  signal.  The  threshold  veloc¬ 
ities  for  spall  formation  are  very  similar  for  1  mm  and  2  mm 
flyers.  For  1  mm  flyers,  plates  uc  =  360  m/s.  For  2  mm  flyer 

plates,  u  =  345  m/s.  Both  of  these  values  are  probably 
c 

accurate  to  ±  10  m/s. 


3.  FREE  SURFACE  VELOCITY  MEASUREMENTS 

Table  12  lists  the  shots  for  which  free  surface  velocity 
measurements  are  available.  Table  13  lists  the  free  surface 
velocity  parameters  for  4340  stael.  The  data  for  the  shots  are 
shown  in  Figures  45,  46,  and  47.  The  oscillographs  for  shot 
219  are  given  in  Appendix  C. 

The  spall  signal  varied  little  with  velocity,  even  though 
in  shot  15  the  release  was  from  the  e-phase.  The  average  value 
of  os  from  equation  (21)  is  -47  kbar. 


4.  SPALL  CRITERIA 


SWAP  calculations  were  carried  out  for  the  threshold 
conditions  of 


(a) 


and 


(b) 


u  =  345  m/s 
c 

d  =2.0  mm 
T  =  4.0  ran 

u  =360  m/s 
c 

d  =  1 . 0  mm 
T  =  5.0  mm 
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A 1  S I  4340  STEEL 
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Figure  44.  Average  Cross  Sectional  Area  of  Voids  Versus 
Strain,  from  COX  and  LOW. 


05  10  1.5  20 

TIME  AFTER  SIGNAL  ARRIVAL  [fis) 


Figure  45.  Free  Surface  Velocity  from  Shot  15. 
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TABLE  12 

FREE  SURFACE  VELOCITY  PARAMETERS  FOR  4340  STEEL 
(Units  =  us  and  mm/us) 


SHOT 


v  . 
nun 


.12  ±  .01 

.15  2  .02 

.52  ±  .01 

.33  ±  .02 

0.77  ±  .02 

.62  ±  .01 

1.14  2  .02 

NR 

1.17  i  .03 

.98  2  .02 

1.65  ±  .03 

1.29  ±  .05 

50  ±  5 

100  ±  10 

175  i  10 

220  +  10 

50  2  5 

52  ±  10 

760  ±  10 

420  2  10 

575  2  5 

200  2  10 

0 .22 

4- 

.02 

0.31 

± 

.02 

0.64 

± 

.01 

0.73 

± 

.01 

1.09 

± 

.01 

1.35 

.L 

.05 

120 

± 

5 

198 

f 

5 

176  , 

106 

2  5 

379 

+ 

5 

180 

+ 

5 

It  was  immediately  apparent  that  these  *  -  v-'rt  incon' 

with  a  Tuler-Bufcchor  approach  because  on  the  ’ '  iane 

o  .  =  -58  kbar  for  care  ia)  but  only  -55  ar  -o  ase  (b)  . 

mm 

(The  a  .  values  in  the  thicker  targets  ai  decreased  by  stress 
min  17 

wave  overtake. 

The  values  of  K~  and  o  used  by  MESCALL  and  PAPIRNO  were 
2  c  o  2 

widely  inappropriate.  For  case  (a)  K-  =  1600  kbar  -us,  and  for 

2  1 

case  (b)  K2  =  500  kbar  -us,  when  cq  -  7  kbar.  For  this  material, 
it  was  decided  to  abandon  the  k,,n  forme t .  The  threshold  data 

A  O 

are  best  described  by  a  single  stress  criterion,  namely 
crg  =  -56.5  kbar.  SWA?  runs  for  shot  219  also  produced  the  best 
agreement  with  Figure  47  for  a  value  of  a  =  ^-60  kbar  (although 
the  uncertainty  is  about  *-i  kbar)-  Thus,  1  340  can  be  described 
by  the  relatively  simple  spall  criterion  cf  og  =  -55  kbar. 
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SEC'1' I  ON  VI 


FRACTURE  OF  ARMOR  PLATE 

The  relied  nomogenous  amor  plat.  v RHA)  tested  in  this 
program  was  RC  33.  It  was  supplied  by  AFATL  in  the  form  of  an 
11  mm  thick  plate.  The  official  specification  for  RHA  are 
given  in  MIL-S-1236QC. 

1 .  PREVIOUS  WORK  ON  RHA 

The  shock  properties  of  RHA  have  been  principally  investi¬ 
gated  by  HAUVEE.  He  found  that  C.  =  5.83  mm/us.  The  HEL  varied 

li 

with  hardness  and  sample  thickness.  His  results  imply  that  for 
the  material  and  specimen  thicknesses  employed  in  this  program, 
the  HEL  should  be  between  22  and  24  kbar.  The  Hugoniot  for 
c.- RHA  was  almcst  the  same  as  a -iron,  namely  U  =  4.51  +  1.43u 
(in  mm/ us ) . 

Fracture  in  RHA  has  been  investigated  by  SHOCKEY  et  al. 

The  NAG  mo *lel  was  used  to  evaluate  flat  plate  and  penetrator 
impacts .  The  material  was  RC38,  which  is  somewhat  harder  than 
that  used  here.  They  found  YQ  =  11  kbar  and  ultimate  strain 
of  0.2  (static  test).  The  stress  threshold  for  dynamic  void 
formation  was  11  kbar. 

2.  SPALL  THRESHOLDS  FOR  RHA 

Table  13  presents  the  spall  threshold  data  for  the  RHA 
tested.  In  all  cases,  separation  was  essentially  complete. 
However,  in  shot  268  the  spall  plate  did  not  completely  detach. 
All  data  were  for  2  mm  flyer  plates.  In  shot  60  there  was  a 
discrepancy  between  the  VISAR  measured  velocity  and  that 
measured  by  the  time  interval  between  crush  pin  closure.  The 
VISAR  data  were  used  in  the  analysis. 
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TABLE  1 : 


SPALL  THRESHOLD  DATA  FOR  RHA 


u 

(m/s) 

d 

(nun) 

T 

(mm) 

shot 

no 

results 

.  . — 

236 

2.0 

4.0 

265 

no  spall 

390 

2.0 

1  4.0 

1 

268 

spall 

771 

2.0 

1 

4.0 

21 

1 

spall 

247 

2.0 

5.0 

62 

no  spall 

380 

2.0 

5.0 

60 

spall 

449 

2.0 

5.0 

123 

spall 

535 

2.0 

5.0 

25 

spall 

622 

_ _ _ i 

2.0 

— 

• 

O 

23 

spall 

For  5  mm  thick  targets,  the  threshold  velocity  is  between 
247  and  390  m/s. 

3.  VISAR  DATA  FOR  RHA 

Good  VISAR  data  were  obtained  in  shots  25,  60,  and  65. 
Limited  VISAR  data  were  obtained  on  shots  21  and  123.  A  sample 
velocity  record,  from  shot  25,  is  shown  in  Figure  48.  Table  14 
presents  the  free  surface  velocity  parameters  determined  from 
the  VISAR  data. 

4  .  SPAT  L  CRITERIA  FOR  RHA 

The  SWAP  calculations  were  used  to  calculate  the  stress 
history  on  the  incipient  spall  planes  for  the  threshold  con¬ 
ditions.  SWAP  calculations  were  carried  out  for  RHA,  using  a 
material  model  consistent  with  the  discussion  in  the  preceeding 


Figure  48.  Free  Surface  Velocity  of  Measured  in  Shot  25. 
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section.  The  threshold  shots  show  that  spall  occurs  at 
og  =  -41  to  -62  kbar  for  d  =  2  mm.  The  spall  stress  calculated 
from  equation  (21)  for  three  shots  that  spalled  (21,  25,  and  60) 
lie  between  46  and  55  kbar,  with  a  mean  of  50  kbar.  It  appears 
a  does  not  change  as  o  varies  from  73  to  150  kbar.  SWAP 

S  IuaX 

calculations  were  not  performed  for  these  shots.  However,  in  a 
series  of  SWAP  calculations  for  spall  in  4340  experiments  with 
this  geometry  (d  =  2  mm,  T  =  4  mm),  a  / ( . 5pCT Au  )  had  a  mean 
value  of  1.74.  Applying  this  result  to  shots  4,  25,  and  60 
gives  a  =  -87  kbar  for  a  wide  range  of  impact  velocities. 

Data  are  not  available  to  evaluate  the  effect  of  stress 
history  on  spall  for  RHA.  The  SWAP  calculations  show  the  stress 
histories  in  4  mm  and  5  mm  thick  targets  struck  by  2  mm  flyer 
plates  are  virtually  identical. 

Based  on  the  data  here,  the  recommended  spall  criteria 

is  a  =  -60  kbar.  This  value  may  slightly  overestimate  threshold 
s 

conditions  and  slightly  underestimate  spall  stress  in  overdriven 
material . 


SECTION  VII 
RESULTS  FOR  NICKEL 


The  material  tested  was  commercially  pure  nickel  200. 

No  previous  spall-related  fracture  data  for  this  material  were 
located.  Hugoniot  data  are  available  from  standard  references, 
namely  VanTHIEL  et  al.,  and  MARSH.  The  low  pressure  Hugoniot 
can  be  represented  as 

U  =  4.59  +  144u 

o 

in  units  of  mm/ys.  The  density  is  aQ  =  8.875  g/cm" . 

CL  =  5.76  mm/ys. 

Spall  threshold  data  for  nickel  are  given  in  Table  15. 

Only  one  VISAR  shot  was  carried  out  on  a  nickel  target,  the 
results  of  which  are  shown  in  Figure  49.  Microstructures  of 
voids  in  nickel  is  shown  in  Figure  50. 

The  VISAR  data  indicate  considerable  residual  stress  on 

the  spall  plane  in  shot  218.  This  is  consistent  with  the 

observation  that  in  229,  which  differed  in  o  by  only  1.5  kbar, 

max 

separation  on  the  spall  plane  was  incomplete.  The  value  of 
a  calculated  from  equation  (21),  using  cn  as  the  spall  shock 
speed,  is  22  kbar.  Using  CL  gives  28  kbar,  but  because  of  the 
low  strength  of  nickel  c0  is  probably  more  appropriate. 

There  is  quite  a  wide  range  of  values  of  impact  velocities 

over  which  uncoalesced  voids  are  observed.  For  2  mm  flyer  plates, 

the  range  is  about  15  Kbar  in  o. _ _  . 

3  max 

There  are  not  enough  data  available  to  adequately  evaluate 
a  time  dependent  spa] 1  threshold  for  nickel.  The  data  available 
imply  a  threshold  stress  for  void  formation  of  20+5  kbar.  This 
applies  for  either  threshold  conditions  or  for  overdriven  spall 
conditions . 
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SECTION  VIII 

RESULTS  FOR  1100  ALUMINUM 


The  aluminum  that  was  the  subject  of  this  study  was 
commercial  grade  1100  alloy.  Sample  microstructure  is  shown 
in  Figure  51. 

1.  PAST  WORK 

The  compressive  properties  of  aluminum  are  well  known, 
and  can  be  found,  for  example,  in  the  compendium  edited  by 
MARSH.  The  Hugcniot  for  1100  aluminum  is  u  =  5.38  +  1.34u 
(mm/ys) . 

The  major  past  study  of  spall  in  1100  series  aluminum  is 
that  of  BARBEE  et  al.  They  found  that  damage  accumulated 
relatively  slowly  over  a  range  of  impact  velocities.  For 
d  =  5.8  mm  and  T  =  15.8  mm,  incipient  damage  occurred  for  flyer 
plate  velocities  of  about  89  m/s.  For  d  =  2.4  mm  and  T  =  6.4  mm, 
the  critical  velocity  was  about  126  m/s.  The  critical  stress 
for  void  nucleation  (in  the  NAG  model)  was  8  kbar.  However, 
the  time  dependence  of  void  growth  and  coalescence  was  quite 
significant. 

2 .  PRESENT  RESULTS 

The  spall  threshold  data  are  shown  in  Table  16.  Complete 
separation  was  very  hard  to  achieve  in  this  material.  Apparently 
the  spall  plate  is  launched  with  such  a  low  velocity  that  shear 
failure  along  its  edges  (necessary  for  complete  detachment) 
could  not  occur.  Figure  52,  for  example,  shows  a  cross  section 
through  the  target  from  shot  51.  Even  though  this  shot  was 
well  above  the  spall  threshold,  a  connected  failure  plane  i^ 
absent.  Spall  plate  thickness  in  aluminum  was  about  25  percent 
less  than  the  flyer  plate  thicknesses.  There  was  a  relatively 
large  velocity  region  in  which  voids  formed  but  a  connected 


Figure  51.  Sample  Preshot  Microstructure  in  Aluminum 


TABLE  16 


SPALL  THRESHOLD  DATA  FOR  HOC  ALUMINUM 


u 

(m/s) 

d 

(nun) 

T 

(mm) 

shot 

no 

results 

300 

1.0 

4.0 

272 

connected  voids 

324 

1.0 

4.0 

271 

partial  crack 

365 

1.0 

4.0 

270 

connected  voids 

472 

1.0 

5.0 

51 

partial  cracks 

590 

1.0 

5.0 

50 

isolated  voids 

1556 

1.0 

5.0 

52 

isolated  voids 

409 

2.0 

5.0 

47 

isolated  voids 

570 

2.0 

5.0 

49 

spall  separation 

894 

2.0 

5.0 

106 

spall  separation 

1024 

2.0 

5.0 

spall  separation 

spall  fracture  system  was  absent.  Figure  53  shows  an  example. 
Spall  occurred  on  three  or  four  parallel  planes,  as  in  the 
model  of  DAVID  et  al. 

SWAP  calculations  were  carried  out  for  the  shots  shown 

in  Table  16.  The  results  for  a  .  are  shown  in  Figure  54. 

mm  3 

The  abscissa  in  Figure  54  is  the  length  of  time  that  the  tensile 

stress  exceeded  n  .  / 2. 

mirr 

The  damage  accumulation  in  aluminum  was  gradual.  The 
spall  criterion  was  taken  as  the  occurrence  of  connected  voids 
as  shown  in  Figure  53. 

The  threshold  for  spall  failure  in  aluminum  was  not  well 

determined  in  these  experiments.  The  impact  velocity  for 

causing  connecting  voids  was  lower  than  could  practicably  be 

obtained  with  the  impact  range  employed.  It  appears  that 

shots  47  and  272  are  close  to  the  threshold.  Calculations 

of  K,  suggest  that  a  value  of  A  =  1  is  much  more  appropriate 

than  \=2.  However,  even  with  X=l,  no  time-dependent  model 

produces  less  spread  between  values  of  shots  27?  and  47 

than  exist  just  in  o  ,  values  for  these  shots.  If  they  are 

both  near  the  threshold,  then  the  best  description  is  simply 

o  <  -20  kbar.  Evidently  over  the  time  scale  of  0.3  to  0.7  ys, 

spall  in  this  aluminum  is  not  time  dependent.  At  a  =  -20  kbar, 

f.  =  -0.028  for  aluminum. 

P 

The  VISAR  data  for  aluminum  shot  272  are  given  in  Figure  55. 
The  spall  signal  is  AV  =  76  ±  5  m/s.  The  residual  stress  on 

9 

the  spall  plane  is  substantial  consistent  with  the  fact  that  the 
voids  did  not  fully  coalesce  in  this  target.  The  value  of  ag 
from  equation  (21)  is  5.5  or  7.4  kbar,  depending  on  whether  the 
bulk  or  elastic  sound  speed  is  employed.  Tnis  is  consistent  with 
literature  values  of  growth  initiation  stress. 
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In  summary,  the  data  for  1100  aluminum  are  not  definitive 
The  most  prudent  interpretation  is  that  significant  void 
growth  begins  in  aluminum  at  valuns  of  os  =  -7.4  kbar,  and  that 
over  time  scales  of  0.3  to  0.7  ys,  spall  formation  is  not 
significantly  time  dependent. 


Figure  55.  Free  Surface  Velocity  from  Shot  272. 
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APPENDIX  A 

EXPERIMENTAL  TECHNIQUES 

The  program  was  carried  out  in  the  Impact  Physics 
Facilities  of  the  University  of  Dayton  Research  Institute. 

These  facilities  are  described  in  detail  by  Hi, ESS  (1981). 

1 •  LAUNCHER 

All  of  the  data  shots  in  the  program  employed  a  40  mm 
solid  propellant  gun.  The  launch  tube  was  4  m  long,  and  is 
composed  of  two  sections.  The  breech  accepts  conventional 
M40  brass  cartridges.  For  the  relatively  low  velocities 
required  in  this  program  an  insert  was  designed  so  that  less 
powder  could  be  used  without  introducing  irregularities  in  the 
loading  curve. 

The  barrel  was  instrumented  with  a  contact  pin  near  the 
muzzle.  A  metallic  ring  on  the  front  of  the  projectile  closed 
the  contact  pin  circuit,  which  generated  a  time  of  arrival 
signal . 

The  projectiles  were  made  from  PMMA.  A  description  of  a 
projectile  is  shown  in  Figure  A-l ,  from  which  most  details  of 
its  construction  are  evident.  For  shots  higher  than  100  m/s, 
the  range  was  evacuated  to  0.5  mm  Hg.  For  lower  velocity  shots, 
the  pressure  was  5  mm  Hg. 

2.  TARGET  MOUNTING 

The  target  was  rigidly  mounted  80  mm  beyond  the  muzzle 
of  the  gun.  Target  orientation  was  adjustable.  Prior  to  each 
shot  perpendicularity  to  the  trajectory  was  verified.  Crush 
pins  provided  time  of  arrival  data  on  the  target.  Initially, 
four  crush  pins  were  used  to  monitor  tilt.  Tilt  was  found  to 
be  always  negligible,  and  in  later  data  shots  only  one  crush 
pin  was  used.  At  velocities  above  200  m/s,  the  pin  was  crushed 
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FLYER  PLATE  a  SABOT 


Figure  a-1. 


Projectile  Design. 


by  the  PMMA  sabot.  At  lower  velocities,  it  was  necessary  to 
strike  the  pins  with  '  ie  flyer  plates  in  order  to  be  certain 
of  prompt  trigger  signals. 

Impact  velocity  was  measured  by  several  different  tech¬ 
niques.  For  shots  with  shot  nuiAber  less  than  212,  impact 
velocity  was  measured  from  two  transit  times.  The  first  was 
between,  a  laser  beam  at  the  muzzle  and  a  crush  pin  12  mm  from 
the  impact  surface.  The  second  interval  was  the  time  between 
closure  of  crush  pins  12  mm  and  1  mm  from  the  impact  surface. 

The  crush  pins  were  manufactured  by  Dynasen  Corporation.  The 
crush  pins  were  positioned  so  that  they  were  struck  by  the  PMMA 
ring  around  the  flyer  plate.  The  average  of  the  velocities  cal¬ 
culated  by  these  techniques  was  employed,  using  inverse  variance 
weighting.  In  212  and  later,  a  Cordin  Model  241  70  mm  streak 
camera  was  used  to  view  the  projectile  pricr  to  impact.  The 
velocity  was  calculated  from  the  photographic  record.  The 
flash  lamp  for  the  streak  camera  was  initiated  from  a  contact 
switch  in  the  barrel  that  was  closed  by  a  conducting  ring 
added  to  the  projectile.  The  velocity  was  also  calculated  from 
the  interval  between  closure  of  the  pin  and  impact  on  a  crush 
pin  1  mm  from  the  target  surface.  High  speed  framing  camera 
pictures  were  taken  during  test  shots  to  verify  the  accuracy 
of  the  velocity  measurement  and  the  integrity  of  the  launch 
package.  For  this  purpose,  a  Beckman  and  Whitley  Model  300 
camera  was  employed.  In  most  shots,  two  or  three  crush  pins 
were  also  dedicated  to  monitoring  impact  planarity.  Figure  A-2 
shows  the  mounting  arrangement. 

The  VISAR  was  mounted  next  to  the  target  tank  on  an  air- 
cushioned  table.  The  VISAR  beam  was  directed  by  three  adjustable 
mirrors.  Initial  alignment  was  done  with  the  target  removed. 

The  VISAR  laser  was  then  adjusted  to  be  coincident  with  the 
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reference  laser  that  shown  along  the  barrel  axis.  At  ATA  VISAR 
was  employed.  The  VISAR  oscilloscope  was  triggered  by  a  crush  pin 
1  or  0.5  mm  from  the  impact  surface.  Sometimes  a  delay  generator 
was  also  used. 


( 2  (LASER  VIEWING 

SLOT  "0"RlNG 


TARGET  MOUNT 


Figure  A- 2.  Target  Mounting  Arrangement. 
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APPENDIX  B 

RELEASE  WAVE  SPEEDS  IN  FERRITE 


Barker  (1974)  concluded  that  in  a  one-dimensional  shock/ 
release  process  in  a  Armco  Iron,  the  release  process  was 
essentially  hydrodynamic.  By  assuming  that  it  is  in  fact 
hydrodynamic,  one  can  calculate  release  wave  speeds. 

An  excellent  description  of  hydrostatic  compression  is 
given  by: 


P  =  3K  f  ( 1- 2 l ) (l+2e£) 
s 


(B-l) 


(see  Birch)  where  K  is  the  isontrop  bulk  modules;  u  is  a  strain 
parameter,  related  to  density  by: 


P  =  PQ  ( l-2c )  '  (B-2) 

and  £  is  a  correction  containing  the  third  order  elastic  content 

C  =  1  «  -  If)  (B-3) 

also 


K  =  K  (1-2:  ) 8/2  [l-7e  +  2£c  (2  +  9e) ] 

b 


( B-4 ) 


where  K  is  here  the  initial  bulk  modules  at  zero  pressure, 
s 

The  bulk  wave  speed  may  be  computed  from: 


CB  =  /k/p  (B- 5 ) 

9  k 

Values  for  these  parameters  are  K  =  1664  kbar,  ttt  a  5.29 

S  dp 

(GUI NAN  and  BESHERS).  Hence  £  =  -0.9675. 

CQ  as  a  function  of  P  can  now  be  computed.  Representative 
values  are: 


P  , 

P 

CB 

c 

(g/cmJ) 

kbar 

km/s 

0.0 

7.85 

0 

4.604 

-0.01 

8.0867 

53.5 

4.726 

-0.02 

8.3257 

114.4 

4.856 
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When  Cn  is  plotted  as  P  between  0  and  130  kbar,  it  is 
seen  tnat  it  is  very  nearly  linear.  The  least-square-best-fit- 
straight  line  is 

Cfi  =  4.6068  +  0.00218P  (B-6) 

This  misses  the  calculated  points  by  only  ~2  m/s  at  both  ends 
of  the  pressure  range. 

If  the  shock  and  release  process  takes  place  wholly  in 

a-iron,  pressure  profiles  may  be  calculated  from  the  method  of 

characteristics.  The  C.  characteristics  are 

+ 

x  =  u  +  C (u) t  (B-7) 


Before  the  release  waves  intercept,  the  C_  invariant  is 


J-  =  uo 

Hence,  anywhere,  and  in  particular  along  C+ 


■dP  _ 
pc 


u  -  u 


(B-8) 

(B-9) 


To  a  very  good  approximation 

/  —  =  -1.25x10"19  P2  +  2.758  x  10"7  P  (B-10) 

J  pc 


when  cgs  units  are  used.  The  relationship  between  A  u  =  u-u0 
and  P  is  then 

P  -  -1103.2  +  /l. 217x10®  -8x10 6 Au  (B-ll) 


where  Au  is  in  km/s  and  P  is  in  kbar  (note  Au<o) . 
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APPENDIX  C 

REDUCTION  OF  VISAR  DATA 


Figure  C-l.  Sample  Visar  Data  Trace 

Referring  to  Figure  Cl? 

A  =  Fringe  amplitude  (contrast) 

Aq  =  Initial  value  of  firing  amplitude 

y  =  Mean  value  of  data  signal 

y  =  Initial  values  of  y 
J  o  2 

h  =  Beam  intensity 

h  =  Initial  value  of  h 
o 

i  =  Subscript  denoting  measurement  time 
The  fringe  amplitude,  f,  is  given  by: 

f  =  (ci) 
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■u,  u*Linw  iinwwflo1 


^w-l1 


or 

0-0 

F  =  - —  (C7) 

360 

If  0  is  expressed  in  degrees.  In  hand-calculator  operations, 
using  0  in  degrees  is  somewhat  simpler  and  the  form  is  used  in 
this  work.  In  BARKER  (1970),  it  is  explained  how  K  is  calculated 
from  the  VISAR  delay  and  laser  wavelength. 

The  most  difficult  stage  of  the  analysis  is  to  determine 
Aq  and  y  .  When  the  peaks  labelled  min  and  max  in  Figure  Cl  are 
discernable,  then: 
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y  -  y  • 

•'max  Jmm 
2 

V  V  . 

'max  +  ■'mm 

2 

Preshot  photos  of  fringe  contrast  were  taken  to  heln  determine 

v  and  y  .  .  However,  it  became  clear  that  during  a  shot  the 

'max  ■ririn  '  ^ 

value  of  Aq  could  be  more  than  in  the  preshot  calibration.  This  is 
probably  due  to  improved  gain  of  the  PM  tubes  on  low  duty  cycles. 

Unfortunately,  much  of  the  time  the  peaks  labelled  min  and 

max  at  the  beginning  of  the  trace  are  not  present,  either  because 

the  shock  rise  time  is  too  fast  for  the  VISAR  to  follow,  or  else 

because  the  maximum  velocity  corresponds  to  less  than  one  full 

fringe.  Usually  late  peaks  are  present.  However,  in  this  case, 

interpretation  must  be  injected  to  distinguish  reversals  from 

peaks.  When  an  unfamiliar  material  is  being  tested,  sometimes 

several  records  must  be  examined  before  reversals  and  peaks  can 

be  definitively  distinguished.  Trial  values  of  A  and  y  must  be 

o  o 

employed  in  order  to  force  f  =  ±1  at  the  designated  points. 

In  most  experiments,  the  peak  free  surface  velocity  is  Known 

from  measurements  of  the  impact  velocity.  This  also  provides  a 

check  on  A  and  y  .  In  some  worst  cases,  it  is  the  only  check  on 
o  o 

these  values;  when  that  happens,  the  VISAR  record  loses  value  as 
a  quantitative  measure  of  velocity,  although  it  still  yields 
accurate  data  for  arrival  times. 

It  is  po  sibie  for  fringe  amplitude  to  change  during  an 
experiment  ir  uependent  of  changes  in  beam  intensity.  This  is 
most  likely  to  occur  when  severe  target  deformation  is  produced 
by  the  shock.  It  is  seldou  a  concern  in  the  present  investigation. 
Changes  in  fringe  contrast  can  be  accounted  for  by  introduction 
of  a  contrast  parameter,  C. 

Initially, 
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At  later  times, 


A. 

1 


C  . 
1 


C.  may  be  still  defined  in  terms  of  C  . 
1  o 


However,  CQ  is  considered  a  slowly-varying  adjustable  parameter. 

Introduction  of  variable  C  can  introduce  considerable  uncer- 

o 

tainty  in  calculated  values  of  u^. 

The  VISAR  records  two  signals  in  quadrative;  that  is, 

separated  in  phase  by  90°.  This  aids  in  the  distinction  of 

fringe  extreme  and  velocity  reversals.  Preshot  calibration 

pictures  were  often  very  useful  in  interpreting  VISAR  records. 

The  preshot  beam  intensity,  h  ,  was  usually  slightly  less  than 

c 

h.  The  value  of  y  could  be  taken  as  equal  to  y  h  /h  .  Often 

o  coo 

A  --Ah  /h  ,  also,  although  sometimes  A  >  A  h  /h  for  reasons 
ococ  ococ 

not  well  understood.  Figure  C-2  shows  a  sample  set  of  data. 
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Figure  C-2 


Two  Quadrature  Channels  from  Shot  219 
(Timing  marks  are  1  MHz.) 
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